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Resumen 



Durante el desarroUo de esta tesis estudiamos el Metodo de Patrones Lummicos Es- 
tandarizables (SCM) para supernovas Tipo II "plateau" haciendo uso de fotometria BVRI 
y espectroscopia optica. Se implemento un procedimiento analitico para ajustar funciones 
a las curvas de luz, de color y de velocidad de expansion. Encontramos que el color V -I 
de estas supernovas, medido hacia el final de la epoca "plateau", puede ser utilizado para 
estimar el enrojecimiento provocado por cl material interestelar de la galaxia anfitriona con 
una precision de (T{Ay) = 0.2 mag. Tras realizar las correcciones necesarias a la fotometria 
se recupera la relacion luminosidad versus velocidad de expansion, reportada previamente en 
la literatura cientifica. Ocupando esta relacion y asumiendo una ley de extincion estandar 
{Rv = 3.1) obtenemos diagramas de Hubble con dispersiones de ~ 0.4 mag en las bandas 
BVI. Por otra parte, si permitimos variaciones en Ry en favor de incertezas menores obten- 
emos una dispersion final de 0.25-0.30 mag, lo que implica que estos objetos pueden entregar 
distancias tan precisas como 12-14%. El valor resultante para Ry es de 1.4 ±0.1, que sugiere 
una ley de extincion no-estandar en nuestra Ifnea de vision hacia este tipo de supernovas. 
Utilizando dos objetos con distancia Cefeida para calibrar la relacion luminosidad- velocidad 
obtenemos una constante de Hubble de 70 ± 8 km s~^ Mpc~^, en buen acuerdo con el valor 
que obtuvo el HST Key Project. 
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Chapter 1 



Introduction 



Supernovae (hereafter SNe) correspond to the explosive, high-energy final stages of some stars. The 
mechanical energy released in these powerful events can reach as much as 10^^ erg (or 1 foe), and their peak 
luminosities can be comparable to the total light of their host-galaxies. SNe can be classified in two types, 
either "Core Collapse" or "Thermonuclear" , depending on their explosion mechanisms. 

Core-collapse SN e (CCSNe) are closely associated to star forming regions in late-type galaxies 



(Anderson fc James 



20081 ). Therefore they have been attributed to massive stars born with > 8 M0 that 



undergo the coll apse of their ir on cores after a few million years of evolution and the subsequent ejection of 
their envelopes (Burrows 200C). Th e se SNe leave a compact object as a remnant, either a neutron star or 



a black hole ( Baade fc Zwickv 



1934 



Arnett 



19961 ). The core-co llapse model receive d considerable support 



with the first detection of neutrinos from the Type II SN 1987A (jSvoboda et al.lll987l ). although no compact 
remnant has been found so far in the explosion site. Among CCSNe we can observationally distinguish 
those with prominent hydrogen lines in their spectra (d ubbed Type II), those with no H b ut strong He lines 



(Type lb), and those lacking H or He lines (Type Ic) (jMinkowski 



1941 



Filippenkolll997l) . Although all of 



these objects are thought to share the same explosion mechanism, their diff'erent observational properties are 
explained in terms of how much of their H-rich and He-rich envelopes were retained prior to explosion. When 
the star explodes with a significant fraction of its initial H-rich envelope, in theory it should display a H-rich 
spectrum and a light curve characterized by a phase of ~ 100 days of nearly con stant luminosity followed by 



a sudden drop of 2-3 mag (jNadvozhin 



2003; 



Utrobin 



2007; 



Bersten et al 



200a). Nearly 50% of aU CCSNe 



belong to this class of Type II "Plateau" SNe (SNe II-P). 

Thermonuclear SNe are characterized by the lack of hydrogen and helium in their spe ctra. Their 
early -time spectra show strong lines due to intermediate mass elements (e.g. Si II, Ca II, Mg II; iFilippenko 



19971) . They are found both in elliptical, spiral, or irregular galaxies. These objects are thought to originate 



in low-mass stars that end their lifes as white dwarfs and explode after a period of mass accretion from a 
companion star, leaving no compact remnants behind them (jHillebrandt fc Niemeveiil2000l ) . Observationally 
they are referred as Type la SNe. 

Given their large intrinsic luminosities, SNe have long been considered potential probes for extra- 
galactic distance determinations and the measurement of the cosmological parameters that drive the Universe 
dynamics. Am ong all types of SNe, the Type la family is the one displaying the highest degree of homogeneity 



( Li et al 



20011) ■ both photometrically and spectroscopically. However, these objects are not perfect standard 
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candles. Empirical calibrations have allowed us to standardize their luminosities to levels of ~ 0.15-0.22 



mag and determine distances to their h ost-galaxies with an unrivaled precision of ~ 7%-10% (jPhillips 



Hamuv et al. 



1996 



1993 



Phillips et al.lll999f l. This powerful technique led a decade ago to the construction of 



Hubble diagrams between z = 0-0.5 and measure very precisely the history of the expansion of the Universe 
over 5 Gyr of look-back time. Contrary to our in tuition these observations revealed that the Universe 



dyna r tiics is described by an a ccelerated expansion (jRiess et al 



1998 



Perlmutter et al 



1999 



Astier et al 



2006 



Wood-Vasev et al.l 120071) . The discovery of the accelerating Universe is profoundly connected with 



theoretical cosmology as it implies the possible existence of a cosmological constant, a concept initially 
introduced by Albert Einstein at the beginning of the 20th century, whose origin still is a mystery. 

Although the acceleration of the Universe has been indirectly c onfirmed by other independent ex- 



perim ents such as the Wilkinson Microwave Anisot ropy Probe (WMAP; Spergel et al 



20031 ) and the Baryon Accoustic Oscillations (BAO; 



Blake fc Glazebrook 



2003 : 



2007 



Bennett et al 



Seo fc Eisenstein 



2003f ). it is 



important to obtain independent con firmation of the Type la results. Although not as bright and uniform as 
the Type la's. lHamuv fc Pintd (|2002l ) showed that the luminosities of SNe II-P can be standardized to levels 
of 0.4 and 0.3 mag in the V- and /-bands respectively, thus converting these objects into potentially useful 
tools to measure cosmological parameters. 

Even though SNe II-P are 1-2 mag fainter and much less homogeneous than SNe la, these objects 
provide two interesting routes to distance determinations. First, the Expanding Photosphere Method (EPM, 
Kirshner fc Kwa nlll974l) . a theoretical technique based on atmosphere models that is independent of the extra- 



galactic distance scale (jEastman et al 



1996 : 



Dessart fc Hillier 



20051 ). This method ca n achieve dispersions 



of 0.3 mag in the Hubble diagr am, which translates into a 14% precision in distance (jSchmidt et al 



1994: 



Hamuv 



2001 



Jones et al. 



2006). Second, the Standa rdized Candle Method (SCM for short), an empirical 



technique initially proposed bv lHamuv fc Pintd (|2002r ) that is based on an observational correlation between 
the absolute magnitude of the SN and the expansion velocity of the photosphere, the Luminosity-Expansion 
Velocity (LEV) relation. This correlation shows that SNe ITP with greater luminosities have higher expansion 
velocities, which permits one to remove the large ('^ 4 mag) luminosity differences display ed by these ob jects 
to levels of only 0.3 mag. So fa r , the SCM has been applied to 24 low-z (z < 0.05) SNe (IHamuvl 120031 ) and 
more recently by 



Nugent et al 



( 20061 ) to 5 high-z (z < 0.29) SNe. The latter work was the first attempt 
to derive cosmological parameters from SNe II-P and demonstrated the enormous potential of this class of 
objects as cosmological probes. 

The upcoming years will witness the deployment of se veral survey telescop es, such as the Panoramic 
Survey Telescope and Rapid Response S ystem fPan-STARRS: iHodapp et al.ll2004l ). the Larg e Synoptic Sur- 
vey T elescope fLSST: Itvsou et al.ll2003l ), the Visible and Infrar ed Survey Telescope (VISTA 



Emerson et al 



2004), the VLT Survey Tele scope (VST; 



Capaccioli et al 



20071 ). and the Skymapper (jCranlund et al 



20031 ). the Dark Energy Survey (DES; 



Castander 



20061) ■ all of which offer the promise to discover SNe by the 



thousands. Whether we use it or not, we will be inevitably confronted by enormous amounts of data on 
SNe II-P which will contain valuable cosmological information. In spite of the great promise shown by the 
SCM, it still suffers from a variety of problems which need to be addressed: 1) the lack of a well-defined 
maximum in the light curves has prevented us from defining the phase of each event; 2) each SN shows 
a different color evolution, which has compromised the use of the photom etric data for th e deter mination 
of host-galaxy extinction; 3) the assumption for dereddening employed bv iHamuv fc Pintd (|2002( ). namely 
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that all SNe reach the same asymptotic color toward the end of the plateau phase, has often led to negative 
extinctions so it still needs to be further tested; 4) the small sample size used so far, especially the scarcity 
of SNe in the Hubble flow, has prevented a proper determination of the intrinsic precision of the method. 

The purpose of this research is to take advantage of the larger and more distant sample of SNe II-P 
available to us today to address the issues mentioned above, refine the SCM, and assess the feasibility of using 
SNe ITP to measure distances in the Universe, in preparation for the massive samples of high-z SNe which 
will be produced in the years to come. With this purpose in mind we have developed a robust mathematical 
procedure to model the light curves, color curves, and velocity curves, in order to obtain a more accurate 
determination of the relevant parameters required by the SCM (magnitudes, colors and ejecta velocities). 
This work makes use of 37 SNe to construct a Hubble diagram (HD), evaluate the accuracy of the SCM, 
and obtain an independent determination of the Hubble con stant. Since our sample has several objects in 



common with the recent EPM analysis of I Jones et al.l (j2008l ). we perform a comparison between SCM and 
EPM. 

We organize this thesis as follows. In § [5] we describe all of the observational material used in this 
work such as the telescopes, intruments and surveys involved. The analysis, methodology and procedures, 
such as the Aq, K and Ahost corrections, are explained in detail in § [3l The dereddening analysis, the 
Hubble diagram, the value of the Hubble constant, and the distance comparison between SCM and EPM are 
addressed in §|31 The final remarks in §[5] provide a discussion about possible variations of the reddening law 
for SNe Il-Palong with comparisons between the value of Hq computed by us and those derived from other 
methods. The first part of this section explores the possibility of a non-standard extinction law in the SN 
host-galaxies. We present our conclusions in § [SI 
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Chapter 2 



Observational Material 



This work makes use of data obtained in the course of four systematic SN follow-up programs carried 
out between 1986-2003: 1) the Cerro Tololo SN program (1986-1996); 2) the Calan/Tololo SN program (CT, 
1990-1993); 3) the Optical and Infrared Supernova Survey (SOIRS, 1999-2000); 4) the Carnegie Type II 
Supernova Program (CATS, 2002-2003). As a result of these efforts photometry and spectroscopy (some IR 
but mostly optical) was obtained for nearly 100 SNe of all types, 51 of which belong to the Type II class 



All of the optical data have been already reduced and they are being prepared for publication (jHamuv et al 



2008[ ). Next we descri be in general t erms the data acquisition and reduction procedures. For more details 



the reader can refer to 



Hamuv et al 



(|2008h . 



2.1 Photometric Data 



The photometry was acquired with telescopes from Cerro Tololo Inter- American Observatory (CTIO), 
Las Campanas Observatory (LCO), the European Southern Observatory (ESO) in La Silla, and Steward Ob- 
servatory (SO). A host of different telescopes and intruments were used to generate this dataset as shown in 
Table 12.11 In all cases we employed CCD detectors and stan dard Johnson-Kron-Cousins-Hamuy UBVRIZ 



filters ( Johnson et al 



1966 : 



Cousins 



1971 



Hamuv et al 



2001) 



The images were processed with IRAI^ through bias subtraction and flatfielding. All of them were 
further processed through the step of galaxy subtraction using template images of the host-galaxies. Photo- 
metric sequences were established aroun d each SN based on observations of Landolt and Hamuy standards 



(Landolt 



1992 : 



Hamuv et al 



1992. 19941 ). The photometry of all SNe was performed differentially with re- 



spect to the local sequence on the galaxy-subtracted images. The transformation of instrumental magnitudes 
to the standard system was done by taking into account a linear color-term and a zero-point. Although this 
procedure partially removes the instrument-to-instrument differences in the SN magnitudes, it should be kept 
in mind that si gnificant systematic discrepancies can still remain owing to the non-stellar nature of the SN 
spectrum (e.g. iHamuv et al.lll990l ). 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities 
for Research in Astyronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Table 2.1 Telescopes and instruments 



Telescope 


Instrument 


x^not/ opec 


\j L i(j u.ym 




T ) 

r 


V4 T O 1 flm 

i.uin 




p 
I 




2DF 


O 


CTIO 1.5m 


CCD 


P 


CTIO 1.5ni 


CSPEC 


S 


Blanco 4.0m 


CSPEC 


S 


Blanco 4.0m 


2DF 


S 


Blanco 4.0m 


CCD 


p 


Swope 1.0m 


CCD 


p 


du Pont 2.5m 


WFCCD 


P/S 


du Pont 2.5m 


MODSPEC 


S 


du Pont 2.5m 


2DF 


S 


du Pont 2.5m 


CCD 


p 


Baadc 6.5m 


LDSS2 


P/S 


Baade 6.5m 


B&C 


s 


Clay 6.5m 


LDSS2 


P/S 


ESO 1.52m 


IDS 


s 


Danish 1.54m 


DFOSC 


P/S 


ESO 2.2m 


EF0SC2 


s 


NTT 3.58m 


EMMI 


s 


ESO 3.6m 


EFOSC 


s 


Kuiper 61" 


CCD 


p 


Bok 90" 


B&C 


s 



Note. — Whether the instrument was used for photometry (P), spectroscopy (S) or both (P/S) is listed in column 3. 
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2.2 Spectroscopic Data 



The spectroscopic data were also obtained with a great variety of instruments and telescopes as 
shown in Table 12.11 The observations consisted of the SN observation immediately follow ed by an arc 



Hamuv et al 



lamp taken at the same position in the sky, and 2-3 flux standards per night from the list of 
|l992. 19941 ). 

We always used CCD detectors, in combination with different gratings/grisms and blocking filters. 
The reductions were performed with IRAF and consisted in bias subtraction, flatfielding, ID spectrum 
extraction and sky subtraction, wavelength and flux calibration. No attempts were done to remove the 
telluric lines. 



2.3 Subsample used for this work 



Of the 51 SNc II observed in the course of these four surveys, a subset of 33 objects comply with 
the requirements of 1) having light curves with good temporal coverage; 2) having sufficient spectroscopic 
temporal coverage; 3) being a member of the plateau class. To this sample we added four SNe from the 
literature: SN 1999gi, SN 2004d.i SN 2004et, and SN 200 5cs. Complementary photometry for SN 2003gd 
obtained bv lVan Dvk et al.l (|2003r ) and lHendrv et al.l (120051 ) was also incorporated in our analysis. Table 
lists our final sample of 37 SNe II-P. For each SN this table includes the name of the host-gal axy, equatorial 



coord inates, the heliocentric redshift and its source, the reddening due to our own Galaxy (jSchlegel et al 
19981 ) and the survey or reference for the data. 
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Table 2.2 Supernova sample 



SN name 


Host Galaxy 


RA(J2000) 


DEG(J2000) 




^host 


(s)b 


E{B- 


-V)gal 


References 


1991al 


LEDA 140858 


19 


42 


24.00 


-55 


06 


23.0 





01525 


HP02 





051 


1 


1992af 


ESQ 340-G038 


20 


30 


40.20 


-42 


18 


35.0 





01847 


NED 





052 


1 


1992am 


MCG-01-04-039 


01 


25 


02.70 


-04 


39 


01.0 





04773 


NED 





049 


1 


1992ba 


NGG 2082 


05 


41 


47.10 


-64 


18 


01.0 





00395 


NED 





058 


1 


1993A 


anonymous 


07 


39 


17.30 


-62 


03 


14.0 





02800 


NED 





173 


1 


1999br 


NGG 4900 


13 


00 


41.80 


+02 


29 


46.0 





00320 


NED 





024 


2 


1999ca 


NGG 3120 


10 


05 


22.90 


-34 


12 


41.0 





00931 


NED 





109 


2 


1999cr 


ESQ 576-G034 


13 


20 


18.30 


-20 


08 


50.0 





02020 


NED 





098 


2 


1999em 


NGG 1637 


04 


41 


27.04 


-02 


51 


45.2 





00267 


NED 





040 


2 


1999gi 


NGG 3184 


10 


18 


17.00 


+41 


25 


28.0 





00198 


NED 





017 


3 


0210 


MGG +00-03-054 


01 


01 


16.80 


-01 


05 


52.0 





05140 


NED 





036 


4 


2002fa 


NEAT J205221.51 


20 


52 


21.80 


+02 


08 


42.0 





06000 


NED 





099 


4 


2002gw 


NGG 922 


02 


25 


02.97 


-24 


47 


50.6 





01028 


NED 





020 


4 


2002hj 


NPMIG +04.0097 


02 


58 


09.30 


+04 


41 


04.0 





02360 


NED 





115 


4 


2002hx 


PGG 23727 


08 


27 


39.43 


-14 


47 


15.7 





03099 


NED 





054 


4 


2003B 


NGG 1097 


02 


46 


13.78 


-30 


13 


45.1 





00424 


NED 





027 


4 


2003E 


MGG-4-12-004 


04 


39 


10.88 


-24 


10 


36.5 





01490 


J08 





048 


4 


2003T 


UGG 4864 


09 


14 


11.06 


+16 


44 


48.0 





02791 


NED 





031 


4 


2003bl 


NGG 5374 


13 


57 


30.65 


+06 


05 


36.4 





01459 


J08 





027 


4 


2003bn 


2MASX J10023529 


10 


02 


35.51 


-21 


10 


54.5 





01277 


NED 





065 


4 


2003ci 


UGG 6212 


11 


10 


23.83 


+04 


49 


35.9 





03037 


NED 





060 


4 


2003cn 


IG 849 


13 


07 


37.05 


-00 


56 


49.9 





01811 


J08 





021 


4 


2003CX 


NEAT J135706.53 


13 


57 


06.46 


-17 


02 


22.6 





03700 


NED 





094 


4 


2003cf 


NGG 4708 


12 


49 


42.25 


-11 


05 


29.5 





01480 


J08 





046 


4 


2003fb 


UGG 11522 


20 


11 


50.33 


+05 


45 


37.6 





01754 


J08 





183 


4 


2003gd 


M74 


01 


36 


42.65 


+15 


44 


20.9 





00219 


NED 





069 


4, 5, 6 


2003hd 


MGG-04-05-010 


01 


49 


46.31 


-21 


54 


37.8 





03950 


NED 





013 


4 


2003hg 


NGG 7771 


23 


51 


24.13 


+20 


06 


38.3 





01427 


NED 





074 


4 


2003hk 


NGG 1085 


02 


46 


25.76 


+03 


36 


32.2 





02265 


NED 





037 


4 


2003hl 


NGG 772 


01 


59 


21.28 


+19 


00 


14.5 





00825 


NED 





073 


4 


2003hn 


NGG 1448 


03 


44 


36.10 


-44 


37 


49.0 





00390 


NED 





014 


4 


2003ho 


ESQ 235-G58 


21 


06 


30.56 


-48 


07 


29.9 





01438 


NED 





039 


4 


2003ip 


UGG 327 


00 


33 


15.40 


+07 


54 


18.0 





01801 


NED 





066 


4 


2003iq 


NGG 772 


01 


59 


19.96 


+18 


59 


42.1 





00825 


NED 





073 


4 


2004dj 


NGG 2403 


07 


37 


17.00 


+65 


35 


58.1 





00044 


NED 





040 


7 


2004ct 


NGG 6946 


20 


35 


25.30 


+60 


07 


18.0 





00016 


NED 





342 


8 


2005CS 


NGG 5194 


13 


29 
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Chapter 3 



Methodology and Procedures 



3.1 AKA corrections 

The photon flux measured by an observer is related to the intrinsic himinosity of the source, its 
distance, dust extinction along the line of sight, and the shift of the spectral energy distribution (SED) to 
longer wavelengths caused by the expansion of the Universe. More specifically, if Ly is the SN rest-frame 
emergent luminosity in units of erg s^^ -^^^ , the photon flux per unit wavelength seen by the observer (in 
units of photons cm^^ s^^ is 



n{X) = ——J (3.1) 

where A' is the SN rest-frame wavelength, A:=A'(1 + z) is the observer's wavelength, and is the luminosity 
distance to the source. 

The flux is modified along its journey to the observer in the following order: host-galaxy reddening 
(Ahost), redshift (i^-term), and Galactic reddening {Aq) (or AKA, for short). In order to be able to extract 
the distance from the observed fiuxes, it proves necessary to remove Nature's imprint on the observed magni- 
tudes in reverse order. To undo Nature's work, we first need to correct the observed magnitudes for Galactic 
extinction (Aq), which is equivalent to moving the observer outside the Milky Way. Then we must move the 
observer just outside the host-galaxy, for which we must correct the spectrum for the redshift caused by the 
expansion of the Universe (K correction) . Finally, we must correct the magnitudes for host-galaxy extinction 
{Afiost), which brings the obs erver to t he SN rest-frame. For the latter step we used as a first approximation 



the reddenings determined bv lDessartI (|2008f ). In a second iteration we applied our own reddenings (see i) l3.3[) . 

In this work the AKA corrections are computed numerically using a library of 196 SNe II-P optical 
spectra. Most of the spectra comes from our own database of 44 Type II SNe, 10 of which are not included 
in Table O (SN 1987A, SN 1988A, SN 1989L, SN 1990E, SN 1990K, SN 1993S, SN 1999eg, SN 2000cb, 
SN 2002gd, and SN 2003ib). These 10 additional objects belong to the Type II class, although do not comply 
with the requirements of § 12.31 to be included in this analysis. The database comprises spectra covering 
the plateau and nebular phases. Each spectrum is brought to the SN rest-frame, i.e., to redshift zero and 
null Galactic reddening using the redshifts and Galactic reddenings listed in Table [^21 Also the spectra are 
corrected using our own host-galaxy reddenings, as determined below (see 13. 3p . 
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3.1.1 Ag corrections 

We define the synthetic apparent Galactic extinction correction as the difference between the mag- 
nitude of an unextinguished SED and the magnitude of the extinguished SED, i.e., 

Ag(X) = -2.5 logio ft 77-^ (3.2 

where 5(A) is tire filter transmission function (see Appendix E]) . We apply this definition to all the spectra 
of our library using the BVRI bandpasses. Since the effective wavelength A changes as the SN spectrum 
evolves, we expect the apparent Aq to be a function of the color of the SED. To examine this point, Figure lXT] 
shows Ag{V) against B — Vioi the specific case = 1 mag, where we identify with different colors 

spectra from the plateau and nebular phases. Although the apparent Aq does not show large variations, the 
dependence on color is evident and must be taken into account. For this purpose wc chose to fit these relations 
with 3rd order polynomials. With this calibration wc can proceed to interpolate the corresponding value of 
the apparent Aq for the specific color of the SN and subtract this value from the observed magnitudes at 
every epoch we have photometry to obtain {mots ^Aq). The RMS of our calibration — 0.0008 mag — provides 
an estimate of the uncertainty in the correction, so we add this number in quadrature to the uncertainty in 
the observed magnitudes. 

3.1.2 K corrections 

We use a similar procedure to compute the K corrections. The synthetic K-term is defined as the 
difference between the magnitude of a redshifted SED and the magnitude of the zero-redshift SED, i.e.. 

Km 9./ /^v A' A^^X') S{X) dX 



This definition is equivalent to that given by ISchneider et al.l (|l983h except that they use fluxes per unit 



frequency. By definition the K-term is a color, therefore we expect it to correlate with the broad-band colors 
of the SED. To demonstrate this point. Figure [32] shows K versus B — V for the specific case of the F-band 
for z=0.05. This dependence is very useful as it permits us to interpolate K values for the specific color of 
the SN and correct the photometry at the plateau or nebular phases separately. The result of this correction 
is the quantity {mobs — Aq ~ K). The RMS of the points around the polynomial fits are of the order of 
~ 0.001-0.01 mag, and these are quadratically added to the error of {viobs — Aq) in order to account for the 
uncertainties in the K correction. 

3.1.3 Ahost corrections 

Finally, we must repeat the previous process for the Ahost correction. The following equation 

A m 9./ !LxXAuost{X) SjX) dX 

A,oM) - -2.5 log.o ji^^^siX) dX ^^'^^ 

defines the apparent Ahost and gives the difference between the magnitude of the supernova right outside the 
host-galaxy {mobs — Aq — K) and the magnitude free of host-galaxy dust {mobs — Aq — K ~ Ahost)- Figure 
shows this correction for all the plateau and nebular spectra of our library separately for the specific case 
^ Aq^'^ is the visual absolute extinction at 5500 A as defined by ICardelli et al.l ll 19891) . 
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Figure 3.1 The dependence of the apparent Ag{V) on the synthetic B — V color obtained from the library 
of spectra must be taken into account. The dashed line defines the input Aq'^'^{V) — 1 mag, indicating 
variations up to 0.025 mag. Plateau-phase spectra (filled circles) are shown and fitted in green, whereas the 
data from the nebular phase are shown with brown triangles. This diagram shows that the Ag{V) versus 
B — V calibration must be treated separately for the plateau and nebular phases. 
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Figure 3.2 The dependence of K{V) on the synthetic color B — V obtained from the hbrary of spectra. 
For a redshift of z = 0.05 the if-terms reach absolute values up to 0.2 mag. Plateau-phase spectra (filled 
circles) are shown and fitted in green, whereas data from the nebular phase are shown with brown triangles. 
The diagram shows that the K{V) versus B — V calibration must be treated separately for the plateau and 
nebular phases. 
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Figure 3.3 The dependence of the apparent Ahost{V) on the synthetic B — V color obtained from the hbrary 
of spectra. The dashed hne defines the input = 1, indicating variations up to 0.023 mag. Plateau- 

phase spectra (filled circles) are shown and fitted in green, whereas data from the nebular phase are shown in 
brown triangles. This diagram shows that the AhostiV) versus B — V calibration must be treated separately 
for the plateau and nebular phases. 
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j4^™j(F]j — 1 mag, which confirms the previous results, namely that the apparent Ahost is a function of the 
broad-band colors. We fit these relations with a 3"^ order polynomial and apply these calibrations to the SN 
magnitudes, making sure to include the RMS of the fits — varying between 0.003-0.01 — in the error budget. 
At this point, we obtain the magnitude of the SN corrected for Galactic extinction, redshift and host-galaxy 
extinction for all epochs (colors) of the SN. 



3.2 Fits to light, color and velocity curves 



In the first incarnation of the SCM 



Hamuv fc Pintd (|2002f) measured all the relevant quantities 



(magnitudes, colors and velocities) at fixed epochs with respect to the time of explosion. In most cases, 
however, it proves hard to constrain this time, thus hampering the task to compare data obtained for different 
SNe. It would be ideal to have a conspicuous feature in the light curves but, unlike other SN types, SNe II-P 
generally do not show an obvious maximum during the plateau phase. 

One way around this is to use the end of the plateau as an estimate of the time origin for each event. 
Although simple in theory, in practice it is not easy to measure this time owing to the coarser sampling of the 
light curves at this phase. Thus, our first aim is to implement a robust light curve fitting procedure in order 
to obtain a reliable time origin to be used as a uniform reference epoch to measure magnitudes, colors, and 
expansion velocities. In the remainder of this section we proceed to implement fitting methods to measure 
reliable colors and expansion velocities. 



3.2.1 BVRI light curve fits 



Figures and [23] show BVRI light curves of eight well-observed SNe II-P. These SNe are repre- 
sentative of the whole sample. As can be seen, there are three distinguishable phases in the light curves: 

— The Plateau phase in which the SN shows an almost constant luminosity during the first ^ 100 days of 
its evolution. This phase corresponds to the optically thick period in which a hydrogen recombination 



wave recedes in mass, gr adually releasing the internal energy of the star (jNadvozhin 



2007 



Bersten et al. 



2008 ) 



2003 



Utrobin 



The Linear or Radioactive Tail, a linear decay in magnitude (exponential in flux) starting about 100 days 
after the explosio n. This phase corresponds to the optically thin period powered by the ^^Co — ^ ^^Fe 
radioactive decay (jWeaver fc Wooslevlll980l ). 



— A Transition phase of ~ 30 days between the plateau and linear phases. 

Both the plateau and linear phases are trivial to model if taken separetely, but the abrupt transition 
makes the fitting task much more challenging, especially with coarsely sampled light c urves. The first a t tempt 



Cleveland et al 



( 19921) 



consisted in using the Local Polynomial Regression Fitting, a technique developed by 
within the framework of the R environment for statistical computing. The name of the method is self- 
explanatory, since it performs a polynomial regression over small local intervals along the domain using a 
routine called loess. Figure [3?6l shows the resulting flts when all the data are fitted simultaneously by loess. 
The small scale features of the plateau are nicely reproduced, but loess is unable to model the transition 



'^^^host visual absolute extinction at 5500 A as defined bv lCardelli et al] ll 19891) . 
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Figure 3.4 Four representative liglit curves of the SN sample. Data points are colored according to the 
filters. BVRI magnitudes are respectively shown in blue, green, brown and red. The dotted line corresponds 
to the analytic fit (see Fig. 13. 8p . In the B and R light curve of SN 2005cs we can appreciate the quality 
of the procedure, since even without data points the code manages to achiev e reasonable fits. The light 
curves of this SN wer e complemented with photometry from the literature (P06: iPastorello et al.ll2006l )(T06: 



Tsvetkov et al 



2006|). 
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Figure 3.5 Four representative light curves of the SN sample. Data points are colored in the same way as 
in Fig. 13.41 The dotted line corresponds to the analytic fit (see Fig. I3.8p . 
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Figure 3.7 Loess fits to the V light curve of SN 1999em (left panel) and to the / light curve of SN 2002gw 
(right panel) . In this case the algorithm was applied separately to each phase (plateau, transition, and tail) . 
Despite the satisfactory accuracy of the regression, we could not find an easy way to force the continuity at 
the two interfaces. 
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satisfactorily. Another attempt was done by fitting separately the three phases of the light curves. The 
results, shown in Figure [3771 satisfy our needs, but the lack of continuity at the two interfaces led us to look 
for alternative approaches. 

After experimenting with several fitting approaches we concluded that the best fits could be achieved 
with analytic functions. After examining several options we ended up using the arithmetic sum of the three 
functions shown in Figure [3751 

> A Fermi-Dirac function (red dashed line in Fig. l3.8p which provides a very good description of the transition 
between the plateau and radioactive phases. 

fFDit) = (3.5) 
1 + e -'0 



ao : represents the height of the step in units of magnitude. 

tpT '■ corresponds to the middle of the transition phase and is a natural candidate to define the origin 
of the time axis. 

wq : quantifies the width of the transition phase. At t = tpT — Bwo the height of the step has been 
reduced by 4.7%, and it decreases down to 95.3% a,t t — tpx + Sivq- 

> A straight line (green dashed line in Fig. 13. 8p which accounts for the slope due to the radiactive decay. 

l{t)^Po{t-tpT) + mo (3.6) 

Po : corresponds to the slope of radioactive tail and an approximate slope for the plateau in units of 
magnitudes per day. 

mo : corresponds to the zero point in magnitude at t = tpT- 

t> A Gaussian function (blue dashed line in Fig. 13. 8p which serves mainly for fitting the B light bump curve 
and the / light curve curvature during the plateau phase. 

g{t) = -Pe-i'^y (3.7) 

The Gaussian function is also useful for reproducing the small scale features that can appear in the 
y-band plateau. 

P : height of the Gaussian peak in units of magnitudes. 
Q : center of the Gaussian function in days. 
R : width of the Gaussian function. 

The resulting analytic function we use to model the light curves, is the sum of the three functions 
detailed above 

nt) = .fFD{t)+m+9{t) 

= +P0 (t - tpT) + mo - Pe<^y (3.8) 

1 + e ™o 
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Figure 3.8 The analytic function used to fit the light curves is shown with the black continuous line. The 
dashed lines represent the decomposition of the main function in its three addends: the Fermi-Dirac in red, 
a straight line in green, and a Gaussian function in blue. 
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which has 8 free parameters. It is fitted to the individual light curves using a ^ minimizing procedure. 
In order to find the minimum we use the Downhill Simplex Method (see Appendix |B| , which although 
is not very efficient in terms of the number of iterations required, provides robust solutions. Examples of 
the analytic fits are shown as dotted lines in Figures 13.41 and 13.51 Although we cannot model most of the 
small scale features in the plateau, the fitting does a really good job modelling the transition. Furthermore, 
the analytic function gives us important parameters that characterize the light curve shape, particularly tpx 
which provides a time origin. In the plots in Fig. 13.41 the time axis is chosen to coincide with the value of 
this parameter obtained from the V light curve. A critical quantity in the analysis that follows is ipT and its 
uncertainty, both of which determine the uncertainties in all the relevant SCM quantities. Normally, when 
the tail phase has been observed the minimizing routine has no difficulties finding tpx and delivers a 
credible error. On the other hand, when the light curve does not have any late-time data points the routine 
underestimates aitpx), in which case we need to provide a more realistic estimate of this uncertainty. The 
criteria to fix (j{tpT) depend on what fraction of the transition phase was sampled. Two useful parameters 
are introduced to describe such sampling: 

tf : the day of the last data point for a given SN. 

tend : the approximate day of the end of the l/-band plateau. We estimate this time around each observed 
point, by calculating the magnitude difference. Am, between the previous and the following point. 
We start from the earliest epochs onward in time, until Am exceeds 0.7 mag. This value is an input 
parameter for the fitting routine (see Appendix [B)) . 

The criteria to estimate tpT and altpx) are the following, 

1. When the transition can be clearly seen, but if we are uncertain that the last data point (at t = tf) 
belongs to the tail phase, we set a conservative estimate of (j(tpT)—5 days. An example of such a case 
is SN 2003hl (Fig.[nH). 

2. When the transition can be clearly seen, but we are sure that the last data point (at t = tf) does not 
belong to the tail phase (see SN 2003hg and SN 2003iq in Fig. 13. 5p , the routine usually fails to converge 
to a reasonable value because there are not enough constraints on the transition phase. We identify the 
failing cases when tpT > tend + 45 because tpT is never greater than tend + 25. In these cases, we set 
tpT — tend + 15 siucc the typical value for the width of the transition is 30 days. For this same reason 
the value of o-(ipT) is set to 15 days, a very conservative value as it encompasses the full possible range 
of tpT values. 

3. When the plateau was extensively observed for more than 100 days but could not detect any hints of 
the transition (see SN 1999br in Fig. 13. 5p . we arbitrarily set tpT = i/ + 15, because the plateau phase 
usually lasts '^100 days. In the two cases where we face this situation, a{tpT) is set to 20 days, which 
generously covers the possibility of a longer plateau. 

Regardless of the sampling of the light curves, we assign a minimum of a{tpT) — 2. 

The light curve fits also have a healthy benefit: they allow us to interpolate magnitudes at epochs 
where only one of the magnitudes was obtained and the second magnitude, necessary to construct a color, 
is missing. Without a color we would be unable to calculate the AKA corrections, so that the interpolation 
feature is extremely beneficial as it permits one to correct magnitudes at all epochs. 
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3.2.2 Color curves fits 

Figure K9\ shows the (B-V), (V-R), and (V-I) colors of three proto-typical Type II-P SNe 
corrected for Aq and i^T-terms. The time origin (the x-axis) corresponds to tpT, i.e. the middle of the 
transition phase. In each case we employ all the data points between day -100 and -10 to fit a Legendre 
polynomial shown with solid lines in Fig. 13.91 The degree of the polynomial was chosen on a case-by-case 
basis and varied between 3rd and 6th order. It is evident that, during the plateau phase, the photosphere 
gets redder with time owing to the decrease of the surface temperature as the SN expands. In theory the 
photospheric temperature should approach and ne ver get below the t empe rature of hydrogen recombination 



around 5,000 K. Based on this physical argument, iHamuv fc Pintd (|2002r ) argued that all SNe II-P should 



reach the same intrinsic colors toward the end of the plateau phase and, therefore, they proposed that the 
color excesses measured at this phase could be at tributed to dust reddening in the SN host-galaxy and be 



exploited to measure Ahost- iHamuv fc Pintd (|2002D performed their analysis with a simple naked-eye estimate 
of the asymptotic colors. Here we improve significantly this situation through the formal color curve fitting 
procedure describe above. 

Armed with the polynomial fits we proceeded to interpolate colors on a continuous one-day spaced 
grid between day -80 and -10 for analyzing colors at multiple epochs (see ? 13. Sp . In a handful of cases the 
data did not encompass the whole grid and we had to extrapolate colors, but never by more than 3 days 
from the nearest data point. 

3.2.3 Fe II based expansion velocity curves 

The third ingredient for SCM is the velocity of the SN ejecta. It is well known that different 
spectroscopic lines yield different expansion velocities. The Fe II A5169 line is thought to closely match the 
SN photospheric velocity and has been usually employed for SCM and EPM. Here we use that line as a proxy 
for the velocity of the SN ejecta. Due to the expansion of the envelope, the spectral lines show a P-Cygni 
profile with an emission centered at the rest-frame wavelength Aq and an absorption shifted bluewards by 
AA. From a measurement of AA we can compute an expansion velocity 

Vexp = C X ^ (3.9) 

Afl 



In this study we use the expansion velocities measured bv I Jones et al.l (|2008l ). Figure 1 01 shows Fe II veloci- 
ties as a function of SN phase, for four different SNe selected for their wide range of sampling characteristics. 
In all cases the SNe show a systematic decrease of their velocities with time. Two physical arguments support 
this observational fact: 1) all the shells of the SN undergo an homologous expansion, i.e. the outer shells 
move faster than the inner shells, and 2) the photosphere recedes in mass allowing us to observe deeper and 
slower layers of the SN as time passes. As shown by Fig. l3.10l the Fe II A5169 expansion velocity curve during 
the plateau phase can be properly modeled with a power law of the form 

Vexp{t) = A X {t - to)" (3.10) 

where A, to and a are three free fitting parameters without obvious physical meaning. 

In general we fit for three free parameters [A, to, a), but when only two velocity measurements are 
available (e.g. SN 1992af in the upper left panel of Fig. l3.10|) we fix the a exponent to -0.5, which corresponds 
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Figure 3.9 {B-V), {V-R), and (V-I) color curves of SN 1999em (filled circles), SN 1991al (open circles) 
and SN 2003hl (filled triangles) corrected for Aq and K-teims. This comparison demonstrates that each SN 
displays a different color evolution, which prevents us to determine color excesses from a simple color offset 
after correcting for Aq and i^T-terms. 
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Figure 3.10 Four expansion velocity curves representative of the SN sample measured by means of the 
Fe II A5169 line profile. The solid lines correspond to a power law (eq. I3.10|) . For SN 1992af (upper left 
panel) we fit only two velocities (see § 13.2.31 for more details). The upper right panel shows a 15 days 
extrapolation for SN 1999cr beyond the last data point. The lower left and right panels show a 10 days 
extrapolation for SN 2003gd and SN 2004dj prior to the first data point. 
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to a typical value for our sample. As shown with solid lines in Figure [3. 101 the fits are quite satisfactory. Here 
we choose to restrict the power-law fits to the plateau phase, since the power-law behavior is not observed 
for expansion velocities beyond the transition phase. We use the same one-day continous grid as the color 
curves in order to interpolate velocities at different epochs between -80 to -10 days. Given the good quality 
of the fits and the shallow slope at late epochs, we allow extrapolations of up to 15 days past the nearest data 
point (see SN 1999cr in Fig. 13.10]) . At the left boundary we reduce the extrapolations to 10 days prior to the 
first point, because the power law gets steeper at early times (see SN 2003gd and SN 2004dj in Fig. I3.10p . 
When only two velocity measurements are available we reduce the extrapolations by 5 days. 



3.3 Host Extinction Determination 

Whil e the d etermination of Galactic reddening is straighforward — thanks to the IR dust maps of 



Schlegeletal 



( 19981 ) — , it is much more challenging to ascertain the extinction due to host-galaxy dust. To 
address this issue here we assume that, owing to the hydrogen recombination nature of their photospheres, all 
SNe II-P should evolve from a hot initial stage to one of similar photospheric temperature. If two SNe have 
similar spectra but suffer different amounts of extinction, all color indices of one object should be redder than 
the corresponding color indices of the other object. This is the hypothesis we want to test in this section. 

Although simple in theory there are a couple of practical difficulties to perform this test. First, it 
is not always possible to contrain the time of explosion and line up color curves from different objects. One 
way around this is to use the transition time t pt defined in § 13.2.11 The second problem is illustrated in 
Figure the color curve shapes can vary significantly from SN to SN, preventing one to measure a single 
color offset between two SNe. Our approach to get around this is to pick a single fiducial epoch in the color 
evolution and assess the performance of such color as reddening estimator. Our polynomial fits to the color 
curves are very convenient for this purpose as they allow us to interpolate reliable colors on a day-to-day 
basis over a wide range of epochs and explore which epoch is the one that gives the best results. 

If all SNe share the same intrinsic temperature at some epoch we expect the subset of dereddened 
SNe to have nearly identical colors (Cq) and the remaining objects should show color excesses, E{C) — C~Cq, 
in direct proportion to their extinctions. A useful diagnostic to check our underlying assumption is the color- 
color plot. Unreddened SNe should occupy a small region in this plane. If we further assume the same 
extinction law in the SN host-galaxies, the subset of extinguished objects should describe a straight line 
originating from such region. The figures of merit in this test are 1) the color dispersion displayed by the 
unreddened SNe, 2) the slope described by the reddened SNe (which is determined by the extinction law), 
and 3) the dispersion relative to the straight line (the smaller the better). 

We have identified four objects in our sample (SN 2003B, SN 2003bl, SN 2003bn, SN 2003cn) 
consistent with zero reddening. Such objects were selected for having: 1) no significant Na l-D interstellar 
lines in their spectra at the redshifts of t heir host - galax ies, and 2) dust-free early-time spectra. For the latter 



we used extinction values determined bv lDessarfi (j2008l ) from fits of Type II-P SN atmosphere models to our 



early-time spectra. Such models use the SN spectral lines to constrain the photospheric temperature and the 
continuum to restrict the amount of extinction. As shown in Figure [3. Ill the atmosphere models successfully 
reproduce the early-time SN spectra. 

We investigated two color-color plots {V-I versus V - R, and V-I versus B — V) over a wide range 
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Figure 3.11 Type II-P SN atmosphere models bv lDessartI (|2008l ) (red line) fitted to our spectra (black line). 
The top panel shows an example of an unsatisfactory fit to one of our late-time spectra. The middle panel 
shows a better fit to an early-time spectrum of SN 1991al. A much better fitting is achieved for an early 
spectrum of SN 2003bl as shown in the bottom panel. 
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Figure 3.12 (V-I) versus (V-R) diagram for 21 SNc II-P having VRI photometry corrected for Aq and 
iiT-terms. The blue hne is a least-squares fit to the data, with a slope of 1.59 ± 0.09. The blue arrow has 
a slope of 2.12 and corresponds to the reddening vector for a standard Galactic extinction law {Ry — 3.1). 
With a blue dot is shown the one SN of this subsample consistent with zero extinction having fi!-photometry. 
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Figure 3.13 {V-I) versus (B — V) diagram for 29 SNe II-P having BVI photometry corrected for Aq and 
iC-terms. The red hne is a least-squares fit to the data, with a slope of 0.77 ±0.04. The red arrow has a slope 
of 1.38 and corresponds to the reddening vector for a standard Galactic extinction law {Rv = 3.1). With red 
dots are shown the four SNe consistent with zero host-galaxy extinction. 
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(from day -50 to -15) of epochs after correcting the photometry for Galactic extinction and X-terms. The 
best results obtained from our scrutiny is the {V - 1) versus (V-R) diagram constructed from day -30 and 
shown in Figure 13.121 At this epoch — approximately the end of the plateau — we obtain the linear behavior 
expected for a sample with the same intrinsic color but different degrees of extinction. Shown with a blue dot 
is the one SN consistent with zero extinction which is, remarkably, one of the bluest objects in this diagram; 
the other three unextinguished SNe do not have i?-photometry. A least-squares fit to the data yields a slope 
of 1.59 ± 0.09, which is close but not exactly equal to the E{V — I)/E(V — R) = 2.12 ratio expected for a 
Galactic extinction law = 3.1), shown as a vector in Figure [3.121 This suggests a somewhat different 
extinction law in the SN hosts compared to the Galaxy. The dispersion of 0.059 in T^-J is a promising result 
as it translates into an uncertainty of AhostiV) — 0.15 mag, which corresponds to the limiting precision 
of this method. The reduced of 1-55 implies that the dispersion can be accounted almost solely by our 
error bars and that any instrinsic color dispersion in our sample is < 0.06. The bottom line is that both the 
(V-R) and (V-I) colors fulfill the minimum requirements as reddening indicators. 

The best results from the (B — V) versus (V-I) analysis were obtained from day -30, which are 
shown in Figure 13.131 The four SNe consistent with zero extinction, shown with red dots, average colors 
{B-V)o = 1.147±0.053 and {V-I)o = 0.656±0.053. Note that there are five SNe in this diagram which are 
slightly bluer than the unreddened sample. The V-I color dispersion of 0.076 is greater than that obtained 
in Fg. 13.121 and is most likely due to the B — V color, since the -B-band is more sensitive to the metallicity of 
the SN, owing to several absorption lines that lie in this spectral region. Therefore we believe that the greater 
dispersion in this diagram could be due to the different metallicities of our SN sample. A least-squares fit to 
the data yields a slope of 0.77 ± 0.04. This slope is quite different than the E{V - T)/E{B - V) = 1.38 ratio 
expected for the Galactic extinction law (shown as a vector in Figure [3.13p . in agreement with the suggestion 
made in the previous paragraph from the {V-I) versus {V-R) diagram. 

We conclude from our exploration that, while the B — V color is problematic, both the {V-R) and 
{y -I) colors offer a promising route for dereddening purposes. In what follows wc will employ solely the 
V-I color since only a small subset of our objects possess R photometry. Although the evidence points to 
a non-Galactic reddening law, for now we will assume a standard reddening law (later on we will relax this 
assumption; see section l4.3.2p . Using our library of SNe II spectra we computed synthetically the appropriate 
conversion factor between E{y — I) and Ay for Type II SNe and a standard reddening law (i?y = 3.1), which 
yielded: 

/3v^^^ = 2.518 (3.11) 

Assuming an intrinsic {V -I)o — 0.656 ± 0.053 the host-galaxy extinction can be computed, with its corre- 
sponding uncertainty, from: 



Av{V-I) = 2.518 X [(y - J) - 0.656] (3.12) 



a{Av) = 2.518 X y'cr(y_7) -I- 0.0532 -I- 0.0592 

where V-I corresponds to the color of a given SN at day -30 (corrected for /-C-terms and foreground 
extinction) and cr(v'-/) combines the instrumental errors in the V and / magnitudes, the RMS of the Aq and 
K terms (see ? 13.1.11 and i; 13.1.21 respectively), and the uncertainty in tpT I3.2.ip . The uncertainty in the 
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Table 3.1 Host-galaxy extinctions for all 37 SNe 



SN name 


^y(spec)^ 


subclass 


Ay(Na I-D)'' 


Aviv-iy 


1991al 


0.31(16) 


silver 


0.31 


(06) 


-0.17(21) 


1992af 


1.24(31) 


coal 


0.17 


;i5) 


-0.37(21) 


1992am 






0.00 


(43) 


0.52(23) 


1992ba 


0.43(16) 


silver 


0.00 


(03) 


0.30(21) 


1993A 


0.00(31) 


bronze 


0.00 


(54) 


0.06(25) 


1999br 


0.25(16) 


silver 


0.00 


(04) 


0.94(25) 


1999ca 


0.12(31) 


coal 


0.34 


(05) 


0.25(21) 


1999cr 


0.47(31) 


coal 


0.69 


(21) 


0.12(21) 


1999em 


0.31(16) 


gold 


1.01 


(05) 


0.24(21) 


1999gi 


0.56(16) 


silver 


0.50 


(08) 


1.02(21) 


0210 


0.31(31) 


bronze 


0.00 


(23) 


0.31(36) 


2002fa 






0.00 


(14) 


-0.35(26) 


2002gw 


0.40(19) 


silver 


0.00 


(02) 


0.18(22) 


2002hj 


0.16(31) 


bronze 


0.00 


(06) 


0.24(22) 


2002hx 


0.16(25) 


coal 


0.00 


(16) 


0.38(22) 


2003B* 


0.00(25) 


silver 


0.12 


(05) 


-0.09(21) 


2003E 


1.09(31) 


coal 


0.71 


(07) 


0.78(23) 


2003T 


0.53(31) 


coal 


0.19 


(19) 


0.35(21) 


2003br 


0.00(16) 


gold 


0.11 


(10) 


0.26(21) 


2003bn* 


0.09(16) 


silver 


0.00 


(03) 


-0.04(21) 


2003ci 


0.43(31) 


coal 


0.00 


(23) 


0.78(23) 


2003cn* 


0.00(25) 


gold 


0.00 


(09) 


-0.04(23) 


2003CX 


0.65(25) 


coal 






-0.27(25) 


2003ef 


1.24(25) 


gold 


1.40 


(12) 


0.98(21) 


2003fb 


0.37(31) 


coal 


0.54 


(23) 


1.24(23) 


2003gd 


0.40(31) 


coal 


0.00 


(04) 


0.33(21) 


2003hd 


0.90(31) 


coal 


0.74 


(27) 


0.01(21) 


2003hg 






2.29 


(12) 


1.97(24) 


2003hk 


0.65(31) 


coal 


1.69 


(20) 


0.44(25) 


2003hl 


1.24(25) 


gold 


1.84 


(09) 


1.72(23) 


2003hn 


0.59(25) 


coal 


0.64 


(08) 


0.46(21) 


2003ho 


1.24(31) 


bronze 


1.28 


(10) 


2.19(21) 


2003ip 


0.40(31) 


coal 


0.42 


(08) 


0.56(22) 


2003iq 


0.37(16) 


silver 


0.91 


(04) 


0.25(22) 


2004dj 


0.50(25) 


coal 


0.26 


(06) 


-0.09(23) 


2004et 


0.00(25) 


coal 


1.17 


(02) 


0.13(27) 


2005CS 










0.72(30) 



Note. — The third column lists the subclass defined upon the criteria exposed in § 14.11 



DessartI 120081) with Ry = 3.1. 



^ equivalent width of Na 1-D measured by us and converted to Ay with a law of 

this research. 
* SNe picked to determine the intrinsic color. 



Barbon et al 



1 19901 1 using Ry = 3.1. 
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intrinsic {V - I)o color was also included in the error of Ay along with the {V -I) RMS in the {V -I) versus 
{V - R) diagram. The host-galaxy reddening values obtained from this technique are listed in column 5 of 
Table \3A \ (with the uncertainties given in parenthesis for the whole sample of 37 SNe). There are eight SNe 
with (V-I) colors bluer than {V -I)o which stand out in this table for their negative reddenings. Although 
not physically meaningful, these negative values are statistically consistent with zero or moderate reddenings. 
In fact, seven out of these eight objects differ by < l.lcr from zero reddening. Only SN 1992af differs by 1.8a 
from Ahost = 0. In the following section we compare this method to other dereddening techniques. 
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Chapter 4 

Analysis 



4.1 Comparing dereddening techniques 



At our request, iDessartI ()2008l ) has kindly performed fits of SNe II-P atmosphere models to our 
library of spectra. In these fits the spectral lines are used to constrain the photospheric temperature and the 
corresponding continuum is employed to estimate the extinction. A crucial condition for this technique to 
work is the spectrophotometric quality of the spectra. 

In general our observations were obtained with the slit oriented along the paralactic angle and the 
relative shape of the spectra should be accurate. However, this was not always possible and sometimes 
the spectra were contaminated from light of the host-galaxy or the slit could not be rotated. For all these 
reasons we first checked the flux calibration of each spectrum by synthesizing magnitudes and comparing 
them to the observed magnitudes, duly interpolated to the time of the spectroscopy. In general we found 
good agreement between the observed and synthetic colors (Figure HH]), thus confirming our confidence in the 
flux calibration. In the 64 cases (18% of all cases) where we found significant differences between synthetic 
and observed colors (> 0.1 for any color) we applied a low-order polynomial correction to the spectrum. 
Basically, this "mangling" correction used the observed photometry to change the slope of a spectrum. After 
checking the flux calibration (and mangling it, if needed) the next step was to correct for Galactic absorption 



and de-redshift the spectra. This database was then used bv iDessarti (|2008l ) for the atmosphere model fits. 
Examples of the fits are shown in Fig. 13.111 

The resulting spectroscopic reddenings are summarized in column 2 of Table [5TT1 As pointed out by 



Dessartj (|2008f) , the spectrum fitting technique works much better with early-time spectra than with late-time 
observations. At late times the photosphere has receded in mass exposing inner and more metal-rich layers, 
which translates into an over-abundance of heavy-elements absorption lines. Therefore the fitting to the 
continuum — practically a temperature fitting — is hampered by the presence of strong absorption lines at 
late times. According to this we divided our sample in four quality subcategories based on the epoch and 
the flux quality of the spectra used in the reddening determination: 

o gold : early-time spectra, without mangling correction 

o silver : early-time spectra, with mangling correction 

o bronze : late-time spectra, without mangling correction 
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Figure 4.1 Synthetic colors minus observed colors versus observed colors for the 150 library spectra. The 
64 cases (18% of all cases) with colors outside the a = 0.1 contours are selected for a "mangling" correction, 
whose purpose is to correct the flux calibration of a spectrum and make it consistent with the observed 
photometry. 
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coal : late-time spectra, with mangling correction 



Note that the main criterion is whether the spectrum is early or late and the second criterion corresponds to 
the flux calibration quality. We trust more the spectra that do not require any corrections as they reflect that 
the observations were better performed, so we consider the unmangled spectra as higher-quality than the 
mangled ones. This sub-classification is given for eac h SN in c olumn 3 of Table [XT] As the reader can notice, 
the error is directly related to this sub-classification. iDessartI (j2008l ) assigns an error of <Je{b-V) ~ 0.05 when 
using early-time spectra, and (JEiB-v) = 0.1 when using late-time spectra. We refine this argument ramping 
up from 0.05 to 0.10, depending on the number of spectra employed for each extinction determination. These 
values are given in Table [5TT] in units of Ay = 3.1 x E{B — V). 

Figure [321 shows a comparison between the spectroscopic reddenings j4v(spec) and our color red- 
denings Av{V — I) derived in ii l3.3[ for the 17 SNe belonging to the top three subclasses (gold+silver+bronze). 
A good agreement is displayed between both techniques with a dispersion of 0.38 mag. The resulting ^ = \A 
suggests that this dispersion is consistent with the combined errors between both techniques. The exceptions 
are two objects: SN 1999br and SN 2003ho. The first object (SN 1999br) only has plateau photometry 
making hard the determination of tpT- Although the error in tpT is quite large (20 days), this uncertainty 
does not have a great impact on the error of the V -1 color due to the flatness of the color curve at these 
epochs (-f 0.0024 per day). Another cause for the disagreement is the extreme properties (low luminosity and 
velocity) of this SN which might also have aV-I color instrinsically redder than that of the bulk of the SNe. 
The second discrepant object (SN 2003ho) belongs to the bronze group so it is possible that the difference 
could be due to the use of a late-time of the spectrum in the determination of the spectroscopic reddening. 

Figure 14.31 shows the same comparison, but this time we include the 16 lowest-quality coal SNe. 
This inclusion clearly deteriorates the good agreement seen in Fig. 14.21 from a — 0.38 mag to cr = 0.51 mag 
(y ^ = 2.0) name ly due to SN 1992af, SN 2003cx, SN 2003hd, and SN 2003fb. This confirms the warning 



bv iDessarti (j2008l ). namely, that his spectroscopic technique works much better with early-time spectra. The 
large value of slso suggests that the errors in the spectroscopic reddenings derived from late-time spectra 
are underestimated. 

Another way to estimate host-galaxy reddening is from the Na 1-D AA5893,5896 interstellar absorp- 
tion doublet observed in the SN spectrum at the host-galaxy redshift. Whenever the line was detected we 
measured its equivalent width; in those cases where we did not detect the Na l-D line we assigned it a null 
value. In all cases we estimated the uncertainty in the equivalent width (EW) based on the signal-to-noise of 
the continuum ar ound this line. We co nverted these measurements into visual extinctions AyiJ^a l-D) using 
the calibration of 



Barbon et al 



mm -. 



E{B ~V) 0.25 X EW{Nal-D), (4.1) 

and tabulate our results in column 4 of Table [3Tl The comparison between Av(Na l-D) and our technique, 
shown in Figure imi exhibits a dispersion of 0.53 mag (x^ = 3.8), which is much higher than the a— 0.38 mag 
obtained from the previous comparison. The disagreement can be attributed to the fact t hat the absorption 
line t races the gas content along the line of sight, but does not necessarily probe dust (jMunari fc Zwitter 
19971 ). Furthermore, reddenings derived from the EW of the Na l-D lines measured from low-dispersion 
spectra simply cannot be expected to be precise which is ex actly our cas e. The basic problem is that the D 
lines produced by a typical interstellar cloud are saturated (jHobbs Ill974l ). The only way that one can hope 
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5 gold 
8 silver 
4 bronze - 



Av (V-1) 



Fig ure 4.2 Comparison between two dereddening techniques: the spectrum fitting method (y-axis, 



Dessart 



20081 ) and our technique (r r-axis, S 13.31). The figure shows three subclasses defined by the quality of the 



spectroscopic data used bv iDessartI (|2008l ): the gold+silver+ bronze sample (see § 14. ip . The solid black line 
is the one-to-one relation. 
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Av (V-I) 



Figu re 4.3 Comparison between two dereddening techniques: the spectrum fitting method (j/-axis, [Dessart 
20081 ) and our technique (x-axis, ij l3.3p . The figure shows the same three subclasses as Fig. l4.2| plus the worst 



quality coal sample (see ii l4.ip . The solid black line is the one-to-one relation. 
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Av (V-I) 



Figure 4.4 Comparison between two dereddening techniques: the Na l-D interstehar hne (y-axis) and our 
technique (x-axis, 13. 3p . The equivalent width of the Na l-D inter stellar line are transformed into visual 

|l99ol) . 



extinctions according to the calibration given by 



Barbon et al 
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to derive the reddening from the D lines is via very high-dispersion spectra that resolve the lines and allow 
the column density to be derived. 



4.2 The Luminosity-Expansion Velocity relation 



Armed with t he dereddening m e thod based on late-time colors we can now revisit the LEV relation 
originally discovered by Hamuv fc Pinto ( 20021) which is at the core of the SCM. For this purpose we applied 
AKA corrections (i) l3.ip to our photometry, we used our analytic fits fi^ l3.2.ip to interpolate BVI magnitudes 
on day -30, and we employed the CMB redshiftj^ in Table WT\ to convert apparent magnitudes to absolute 
values (assuming Hq = 70 km s^^ Mpc~^). The expansion velocities were determined from the minimum of 
the Fe II A5169 P-Cygni line profiles. We performed a power-law fit to interpolate a velocity contemporaneous 
to the photometry (day -30) as described in section 13.2.11 From our original sample of 37 SNe II-P we were 
able to use 30 SNe to build this relation, since five of them do not have Fe II velocities at day -30, and two 
others have extremely low redshifts {czcmb < 300 km s^^). 

Figure [4751 shows the resulting LEV r elation (absolute magni tude versus expansion velocity) for all 
BVI bands. Evidently we recover the result of lHamuv fc Pintd (j2002r ). namely that the most luminous SNe 
have greater expansion velocities. In their case the data were modeled with a linear function. Our sample 
suggests that the relation may be quadratic, but we need more SNe at low expansion velocities to confirm 
this suspicion. Linear least-squares fits to our BVI data yields the following solutions 



Mabs{B) = 3.50(±0.30) log(uFcii/5000) - 16.01(±0.20) 
MabsiV) = 3.08(±0.25) log(uFcii/5000) - 17.06(±0.14) 
Mabsil) = 2.62(±0.21) log(uFcii/5000) - 17.61(±0.10) 



(4.2) 
(4.3) 
(4.4) 



which are shown with solid lines in Figure 14.51 The relation found by iHamuv fc Pintol (|2002l) for the y-band 



i s show n with the dashed line in the middle panel of the same figure. Given that the study of lHamuv fc Pinto 
(|2002l ) was performed using data at day 50 after the explosion (around day -60 in our time scale), it is not 
unexpected that their LEV relation is shifted to higher expansion velocities. Some of the difference in slope 
is explained by the inclusion of SN 2003bl in our sample, which flattens the correlation. This relation exhibits 
a dispersion of 0.3 mag, similar to that reported before. This low scatter is very encouraging as it implies 
that the expansion velocities can be used to predict the SN luminosities, to standardize them, and to derive 
distances. 



^the CMB redshifts were computed by adding the heliocentric redshifts given in Tabic [2?2l and the projection of the velocity 
of the Sun relative to the CMB in the directi on of the SN host galaxy. For the latter wc adopted a velocity of 371 km s^^ in 
the direction {I, b) = (264.1, 48.3) given bv iFixsen et al.l lIlQQd) . 
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Figure 4.5 BVI band absolute magnitudes (j/-axis) against the Fe II expansion velocity (x-axis) for 29- 
30 SNe. We use the magnitudes corrected for Galactic extinction, A'-terms, and host-galaxy extinction 
together with a value of the Hubble constant of 70 km s~^ Mpc~^. The d ashed line in the middle panel 
represents the LEV relation for the V^-band found bv lHamuv fc Pintd (|2002r ) obtained from magnitudes and 
velocities measured at day 50 past the explosion (approximately day -60 in our own time scale). 
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4.3 Hubble diagrams 



The LEV relation shown in the previous section imphes that a spectroscopic measurement of the 
expansion velocity of a SN II-P can be used to compute a corrected luminosity which should be approximately 
the same for all SNe. We examine the reality of this property of SNe II-P in the magnitude-redshift Hubble 
diagram. For this purpose we employ Fe II velocities (in units of km s^^), BVI apparent magnitudes corrected 
for K-teims, Galactic reddening, and host-galaxy reddening determined from V-I colors, and host-galaxy 
redshifts in the CMB frame. Table |4?T] lists these values for the 37 SNe of our sample, of which 35 meet the 
requirement of being in the Hubble flow {czcmb > 300 km s^^). A perfect distance indicator would describe 
a straight line in the Hubble diagram, so the figure of merit to assess the performance of this method is the 
dispersion from the fit. Along this section we use dispersions weighted by the errors to evaluate the precision 
of the Hubble diagrams. 



4.3.1 Using Av{V-I) and Ay(spec) 

The top left panel of Figure 14.61 shows a Hubble diagram constructed from _B-band magnitudes 
interpolated to day -30 previously corrected for Galactic extinctions and /C-terms. The top right panel shows 
magnitudes additionally corrected for the LEV relation, the bottom left panel includes further corrections 
for host-galaxy extinction (using the V - 1 color calibration given in section 13. 3p . In each case we perform a 
linear least-squares fit of the form. 



m + a log(i7Fcii/5000) - Ahost = 5 log cz + zp (4.5) 

where m is the apparent magnitude corrected for if-tcrms and Galactic absorption, vpeii is the expansion 
velocity, Ahost is the host-galaxy absorption, and z is the CMB redshift. The only fitting parameters are a 
and zp; in the top left panel we set a = and we only fit for the zero point. 

The dramatic decrease in the dispersion, from 0.72 to 0.36 mag clearly demonstrates the benefitial 
effects of adding the velocity and host-galaxy extinction terms. An inspection of the T/-band diagram fFig. l4.7p 
shows a large scatter of 0.54 mag in the top left panel. When we correct for expansion velocities, the scatter 
drops to 0.50 mag. This is certainly not unexpected given the LEV relation reported in the previous section. 
It is encouraging to notice that the dispersion drops from 0.50 to 0.45 mag when we include our host-galaxy 
extinction corrections. This indicates the usefulness of our color-based dereddening technique. The reduced 
value of 2.45 implies that most of the scatter is accounted by the observational errors. We performed the 
same analysis using other epochs and we found that day -30 yielded the lowest dispersion. This is also the 
day for which we reach the maximum number of SNe in our HDs, i.e. moving backwards or forwards in time 
means losing SNe data (velocity or magnitudes out of the observation range). 

If we turn our attention to the /-band (Figure 14. 8p the final dispersion is 0.45 mag, identical to 
that found in the y-band. It seems that the dispersion could have some dependence on wavelength, since 
it decreases from 0.45 in VI to 0.36 mag in B. However, it may be due to a sampling effect, because the 
VI diagrams have one SN more than the B diagram. The scatter of ~ 0.4 mag in BVI is comparable but 
somewhat larger than the ^ 0.34 dispersion obtained in previous SCM studies inamuv fc Pint J2002[ 



Hamuv 



20031 ). It is important to notice that the dispersions are independent of the intrinsic {V-I)o color calculated 



in section 
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Table 4.1 Magnitudes, expansion velocities, and V -I colors for day -30 



SN name 


CZCMB* 




B 






V 






J 




WFoII 




V 


/ 


1991al 


4480 


17 


944 


^036) 


16 


917 


;029) 


16 


333 


^020) 


5328(202) 





589 


^022) 


1992af 


5359 


18 


024 


;038) 


17 


096 


;029) 


16 


706 


;053) 


4529(212) 





511 


;027) 


1992am 


14007 


20 


443 


'117) 


19 


053 


;038) 


18 


218 


;033) 







830 


'046) 


1992ba 


1245 


16 


962 


^040) 


15 


657 


;o3i) 


14 


886 


^022) 


2237(068) 





776 


^021) 


1993A 


8907 


20 


636 


;i27) 


19 


309 


;089) 


18 


657 


;057) 







678 


;047) 


1999br 


1285 


19 


010 


;022) 


17 


582 


;oi3) 


16 


574 


;008) 


1127(205) 


1 


028 


'Oil) 


1999ca 


3108 


17 


901 


;o7i) 


16 


489 


;055) 


15 


735 


;034) 







755 


'025) 


1999cr 


6363 


19 


435 


;067) 


18 


418 


;o5i) 


17 


723 


;035) 


3095(207) 





703 


;025) 


1999cm 


670 


15 


331 


;034) 


13 


998 


;023) 


13 


245 


'017) 


2727(056) 





752 


;oi9) 


1999gi 


831 


16 


554 


;032) 


15 


060 


:oi9) 


14 


Oil 


;022) 


2725(061) 


1 


058 


;025) 


0210 


15082 


21 


955 


;044) 


20 


572 


;028) 


19 


733 


;035) 


4923(206) 





780 


;iio) 


2002fa 


17847 


21 


222 


;072) 


20 


243 


:052) 


19 


709 


;088) 


4176(061) 





518 


;065) 


2002gw 


2878 


18 


575 


;020) 


17 


491 


:029) 


16 


752 


;oi4) 


2669(063) 





726 


;035) 


2002hj 


6869 


19 


705 


;084) 


18 


582 


:o6i) 


17 


939 


;042) 


3657(065) 





751 


;031) 


2002hx 


9573 


20 


328 


;056) 


19 


164 


;035) 


18 


360 


;028) 


3960(203) 





805 


;o3i) 


2003B 


1105 


17 


078 


;022) 


15 


966 


;oi8) 


15 


341 


;oi5) 


2795(067) 





620 


;oi7) 


2003E 


4380 


19 


659 


;046) 


18 


368 


;026) 


17 


472 


^024) 


2869(203) 





965 


;024) 


2003T 


8662 


20 


541 


;028) 


19 


227 


;o2o) 


18 


423 


^022) 


2803(203) 





795 


^024) 


2003bl 


4652 


20 


158 


;042) 


18 


961 


;o2o) 


18 


230 


'019) 


1475(202) 





758 


'025) 


2003bn 


4173 


18 


638 


;045) 


17 


403 


;034) 


16 


769 


;024) 


2719(202) 





641 


;020) 


2003ci 


9468 








19 


634 


;035) 


18 


648 


^029) 


2876(077) 





963 


'034) 


2003cn 


5753 


19 


852 


;024) 


18 


827 


;oi5) 


18 


182 


;oi4) 


2510(207) 





640 


'017) 


2003CX 


11282 


20 


417 


;ioo) 


19 


513 


;o5o) 


19 


050 


;067) 


3734(205) 





551 


;058) 


2003ef 


4504 


19 


304 


;032) 


17 


837 


;024) 


16 


792 


;017) 




1 


043 


;oi5) 


2003fb 


4996 


20 


521 


'122) 


19 


085 


;092) 


17 


940 


;056) 


3120(206) 


1 


148 


;040) 


2003gd 


359 


15 


208 


'045) 


13 


965 


:036) 


13 


167 


;023) 


2976(230) 





786 


'019) 


2003hd 


11595 


20 


462 


;031) 


19 


320 


:014) 


18 


658 


'018) 


3741(069) 





661 


;022) 


2003hg 


3921 


20 


247 


;091) 


18 


067 


;o4o) 


16 


603 


'025) 


3398(211) 


1 


435 


;027) 


2003hk 


6568 


19 


478 


;036) 


18 


201 


;024) 


17 


380 


;023) 


3618(087) 





828 


;026) 


2003hl 


2198 


19 


124 


;062) 


17 


219 


;038) 


15 


806 


;027) 


2354(063) 


1 


336 


;027) 


2003hn 


1102 


16 


416 


;020) 


15 


153 


;oii) 


14 


307 


'015) 


3121(074) 





836 


'018) 


2003ho 


4134 


20 


779 


;084) 


18 


946 


;026) 


17 


419 


;016) 


4152(080) 


1 


523 


;022) 


2003ip 


5050 


18 


894 


;046) 


17 


549 


;035) 


16 


664 


;025) 


3813(204) 





876 


'021) 


2003iq 


2198 


17 


396 


;049) 


16 


124 


;037) 


15 


362 


;024) 


3482(063) 





756 


'018) 


2004dj 


180 


12 


973 


;074) 


12 


046 


;035) 


11 


416 


;033) 


2725(202) 





621 


;043) 


2004et 


-133 


13 


683 


;220) 


12 


099 


;i7i) 


11 


378 


^104) 


2901(202) 





707 


^070) 


2005CS 


635 


16 


099 


;062) 


14 


749 


;o5i) 


13 


824 


^055) 







957 


^067) 



* The Velocity Calculator tool at the NASA/IPAC Extragalactic Database webpage computes the CMB rcdshift (czqmb in 
km s~^) given the coordinates and the hehocentric redshift of an object. The error of this quantity is dominated by the error 
in the determination of the Local Group velocity, which is 187 km s~^. 

Note. — All the data in this table has been interpolated to day -30 with errors accounting for the interpolation, the error 
oitpT, and the instrumental error. BVI corresponds to apparent magnitudes corrected for Galactic absorption, if -terms, but 
uncorrected for host-galaxy extinction. The V-I colors do not exactly match the V minus / magnitude difference, since light 
and color curves are interpolated independently. 
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Figure 4.6 S-band Hubble diagrams. The top panel shows magnitudes interpolated to day —30 previously 
corrected for Galactic extinction and iiT-terms; the top right panel shows magnitude additionally corrected for 
expansion velocities; the bottom left panel includes further corrections for color-based host-galaxy extinction 
(using the V-I color calibration given in s ection 13.31): th e bottom right panel replaces the color-based 
extinction for the spectroscopic reddenings of Dessart pOOS ) . 
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Figure 4.7 y-band Hubble diagrams. The top panel shows magnitudes interpolated to day -30 previously 
corrected for Galactic extinction and /-C-tcrms; the top right panel shows magnitude additionally corrected for 
expansion velocities; the bottom left panel includes further corrections for color-based host-galaxy extinction 
(using the V-I color calibration given in s ection 13.31): th e bottom right panel replaces the color-based 
extinction for the spectroscopic reddenings of Dessart pOOS ) . 
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Figure 4.8 /-band Hubble diagrams. The top panel shows magnitudes interpolated to day —30 previously 
corrected for Galactic extinction and iiT-terms; the top right panel shows magnitude additionally corrected for 
expansion velocities; the bottom left panel includes further corrections for color-based host-galaxy extinction 
(using the V-I color calibration given in s ection 13.31): th e bottom right panel replaces the color-based 
extinction for the spectroscopic reddenings of Dessart pOOS ) . 
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In the bottom rig ht panel s of Fi gs. 14.61 14.71 and 14.81 we examine the SCM using the spectroscopic 
reddenings determined by iDessartI pOOST ) for a set of 28 SNe II-P (Table 13. ip instead of our color-based 
extinctions. The resulting dispersions in BVI are (0.67, 0.54, 0.42), which compare to (0.36, 0.45, 0.45) 
when using color-based extinctions. All the fitting parameters derived from both dereddening techniques are 
compiled in Table [ 



4.3.2 Leaving Ry as a free parameter 

The previous analysis suggests that the dispersions are somewhat larger than the observational 
errors. One possible source of scatter could be the reddening law. As shown in 13.31 and Fig. 13.121 we have 
evidence pointing toward a somewhat different extinction law in the SN hosts compared to the Galaxy. Here 
we take this idea a step further and we analyze the Hubble diagram leaving Ry as a free parameter. To 
accomplish this we model the data with the following expression 



m-Ka log(uFcii/5000) -/3(F-/) = 5\ogcz + zp (4.6) 



Note that we are replacing Ahost hi equation 14.51 with the term f3{V-I). The /3 factor is a new free pa- 
rameter to be marginalized along with a and zp, and V - 1 is the color on day -30 corrected for _fC-terms, 
Galactic extinction but uncorrected for host-galaxy dust. Once /3 is known we can used it to solve for host- 
galaxy extinction from A{X) = (3 (V^-/)— 0.656 , where 0.656 is the intrinsic V -I color of SNe H-P (see 
equations 13.111 and 13.121 in § 13. 3|) . 

Figure 14.91 shows the BVI Hubble diagrams for our set of 30 SNe H-P. We get dispersions of 
(0.28, 0.31, 0.32) in BVI respectively, which compare to (0.36, 0.45, 0.45) when Ry is kept fixed. The 
increase in is due to the fact that we are not using the intrinsic color, which gives the major contribution 
to the errors. Although we expect a reduction in the scatter due to the inclusion of an additional parameter, 
the large drop in the dispersion is remarkable. The last three lines of Table 14.21 shows the parameters we 
obtain by minimizing the dispersion in the HD for each band. When we restrict the sample to objects with 
czcMB > 3000 km s~-^ (leaving aside SNe with potentially greater peculiar velocities), we end up with 19 SNe 
in the S-band and 20 SNe in the VI bands. The resulting HDs in BVI show dispersions of (0.25, 0.28, 0.30), 
respectively. Not surprisingly these are lower than the (0.28, 0.31, 0.32) dispersions obtained from the whole 
sample, and the Hubble constant shows only a mild increase of 3%. 

By definition the term [3{V-I) in equation 14.61 corresponds to the extinction in a broad-band 
magnitude (with central wavelenth A). Thus, j3 is the ratio A{X) / E{V ~ I) (see eq. l3.11l in ii l3.3p and is related 
to the shape of the extinction law. For each of the BVI bands we used our library of SNe H spectra to compute 
synthetically the value of j3\ as a function of Ry ~ Ay /E{B — V) (see Figure l?.10l) . This allowed us to convert 
the /3a values resulting from our fits into Ry values. Our fits yield j3x = (2.67±0.13, 1.67± 0.10, 0.60 ± 0.09) 
for the BVI bands, which translate into Ry — (1.38lg'24, l-44io'i4, l-36^g :[2)- These values are remarkably 
consistent to each other and significantly lower than the Ry = 3.1 value of the standard reddening law. 
Independent evidence for a low Ry law has been already reported from studies of SNe la (see § [5]) . 

The dispersion of 0.28-0.32 mag in the HDs translates into an accurac y of 13-14% in the determi 



nation of distances. Altho ugh not as good as the ^ 7%-10% precision of SNe la (jPhillips 



1996 : 



1993 : 



Hamuv et al 



Phillips et al.l 119991 ). this is a very encouraging result which demonstrates that SNe II-P have great 



potential to determine extra-galactic distances, and therefore, cosmological parameters. 
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Figure 4.9 BVI Hubble diagrams leaving Ry as a free parameter. 
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Table 4.2 Fitting parameters from the Hubble diagrams 



Reddening 










estimator 


filter 


a 




zp 


(V-I)^ 


B 


3.27 ±0.37 




-0.24 ±0.08 




V 


3.01 ±0.36 




-1.36 ±0.08 




I 


2.99 ±0.36 




-2.02 ±0.08 


spectral 


B 


4.77 ±0.47 




-0.60 ±0.11 


analisis^ 


V 


4.22 ±0.45 




-1.60 ±0.10 




I 


3.67 ±0.43 




-2.11 ±0.10 


(V-I) + 


B 


3.50 ±0.30 


2.67±0.13 


-1.99 ±0.11 


variable Ry'^ 


V 


3.08 ±0.25 


1.67±0.10 


-2.38 ±0.09 




I 


2.62 ±0.21 


0.60 ±0.09 


-2.23 ±0.07 



" Reddenings from § 13.31 
^ Reddenings from DessartI ||2008j) 



Leaving Ry as a free parameter 
Note. — These are the results of fitting 

m + a log (?;pgji/5000) - A^ost = 5 log cz + zp 

and 



m + a log (?jp(,ii/5000) - I3(V-I) = 5 logcz + zp 



to the SN data in Tables [Q and Ol 
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Figure 4.10 Ry versus f3. The solid lines are computed from our library of SNe II spectra for each of the 
BVI bands. The f3 parameters derived from the HDs in Figure [4?9l are shown with red open circles for each 
of the BVI bands. The red dotted line shows the Ry weighted mean for the BVI values of f3. The values 
for the standard Galactic extinction law are also shown with green dots. 



46 



4.4 The Hubble constant 



The Hubble constant is a parameter of central importance in cosmology which can be determined 
from our Hubble diagrams. This can be accomplished as long as we can convert apparent magnitudes into 
distances. This calibration wa s done with two obje cts for which we found Cepheid dis tances in the literature : 
SN 1999em {ji = 30.34 ±0.19; 



Leonard et al. 



20021 ) and SN 2004dj {^l = 27.48 ± 0.24; 



Freedman et al 



20011) 



For each calibrating SN we can solve for Hq using the following expression 



log iJo = 0.2 X m + a log (ufoIi/5000) - f3{V - 1) - n - zp + 25 



(4.7) 



In this equation m is the apparent magnitude corrected for Galactic absorption and iiT-terms, vpeii is the 
expansion velocity, and - / is the color of the calibrating SN, all measured on day -30 and given in Table HTTl 
a, /3, and zp are the fitting parameters given in the bottom three lines of Table 

Table 14.31 summarizes our calculations for the two calibrators. This table shows that the corrected 
BVI absolute magnitudes of SN 1999em and SN 2004dj differ from each other by 0.86-1.13 mag. This is 
somewhat greater than the scatter of 0.3 mag in the LEV relation, but statistically plausible. This can be 
seen in Figure 14.111 where we plot with red dots the corrected absolute magnitudes of the two calibrating 
SNe, on top of the whole sample of SNe employed in the HDs (black dots). This difference leads to Hq values 
in the range 62-105 km s~^ Mpc~^. A weighted average of the values of Table vields a Hubble constant 



i7o = 70 ± 8 km s^^ Mpc'^ 



(4.8) 



This value compares very well with tha t derived by 
and reasonable well with that found by 



Freedman et al 



Sandage et al 



pOOll ) from SNe la (71±2 km s'^ Mpc'^), 



20061 ) (62.3±1.3 km s"^ Mpc"^) with a similar Type la 



sample. With only two calibrating SNe the SCM still has plenty of room to deliver a more precise value for 
Ho. 



Table 4.3 Hq calculations 





M + a 


log(u/5000) - f3{V - I) 








SN 


B 


V I 


Ho{B) 


Ho{V) 


Hoil) 


1999em 


-17.94±0.20 


-18.41 ±0.20 -18.24±0.19 


64.7 ±6.8 


62.2±6.1 


62.9 ±6.0 


2004dj 


-17.08 ±0.30 


-17.28 ±0.27 -17.13 ±0.26 


95. 7± 14.0 


104.5 ±13.7 


104.8 ±12.9 


Average 


-17.67±0.39 


-18.02±0.53 -17.84±0.53 


71 ±12 


69 ±16 


70 ±16 



Note. — The uncertainties in (Hq) and {Mcorr) (the corrected absolute magnitude for BVI in the second column) correspond 
to the dispersion in each pair of values, which are about x2-3 greater than the formal errors. 
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Figure 4.11 BVI corrected absolute magnitudes (y-axis) against the Fe II expansion velocity (a;-axis) for 
29-30 SNe (black dots). We use the Ho values of Table and the CMB redshifts to compute the distance 
moduli. The dashed lines show the mean corrected absolute magnitude for the black dots. The red dots are 
the two calibrating SNe, whose corrected absolute magnitudes were calculated using their Cepheid distances 
(see Table . Note that, as expected, the two calibrating SNe fall on each side of the corrected absolute 
magnitude distributions. 
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4.5 Distances 



Armed with standardized absolute magnitudes for SNe II-P, we are in a position to calculate dis- 
tances to all the SNe of our sample. This can be done with the following expression 



/i = rn + a log (ufcIi/5000) - f3{V-I) - {Mcorr) 



(4.9) 



where {Mcorr) is the weighted mean of the corrected absolute magnitudes of the two calibrating SNe given 
in Table 14.31 The resulting values are given in Table 14.41 for the 37 SNe II-P. The last column of the table 
shows the weighted averages of their distance moduli. 

As mentioned in the introduction, we can evaluate the precision of the SCM from a comparison 
with EPM distances. For this purpose we employ EPM distances recently calculated by I Jones et al.l (j2008r ) 
which are summar i zed in Table B31 along with our distances for the 11 objects in common between SCM and 
EPM. 



Jones et al 



( 20081) determine EPM distances with two di fferent sets of atmosphere models. The EPM 



distances determined using the models of 



Eastman et al. 



( 19961) (E96; column 2) are 12% d z 5% shorter than 



Dessart fc Hilliei 



our SC M distances. On the other hand, the EPM distances determined using the models of 
( 20051) (DOS; column 3) are 40% ± 10% greater than our SCM distances. These shifts are calculated weighting 
by the error in the EPM and SCM distances. These systematic differences can be clearly seen in the upper 
panel of Figure [4.121 and more clearly in the fractional differences ((Iepai — dscAi) / dscM plotted in the 
bottom panel. 

The systematic differences among the two sets of EPM distances can be solely attributed to the 
atmosphere models of E96 and D05. New radiative transfer models of SNe II-P are necessary to understand 
the origin of this discrepancy. Besides the systematic errors in both EPM implementations, we can get 
an understanding of the internal precision of EPM and SCM after removing the systematic differences and 
bringing all distances to a common scale. For this purpose we correct the EPM distances to the SCM 
distance scale by removing the percentage shifts of 1.40 and 0.88. The comparison is shown in Figure [4. 131 
The distance differences have dispersions of ~ 13% and ^ 16% using DOS and E96 respectively. This implies 
that either SCM or EPM produce relative distances with an internal precision between 13-16%. This agrees 
with the dispersions of 0.2S-0.30 mag seen in the HDs restricted to SNe with czcmb > 300 km s~^. 
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Table 4.4 Distance Moduli 
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2003gd 


M74 


29 


99(42) 


29 


98(55) 


29 


94(54) 


29 


98(28) 


2003hd 


MCG-04-05-010 


35 


93(40) 


35 


85(54) 


35 


77(53) 


35 


86(28) 


2003hg 


NGC 7771 


33 


50(42) 


33 


18(54) 


33 


14(54) 


33 


31(28) 


2003hk 


NGC 1085 


34 


45(40) 


34 


41(54) 


34 


35(53) 


34 


41(28) 


2003hl 


NGC 772 


32 


08(41) 


32 


01(54) 


31 


10(53) 


32 


04(28) 


2003hn 


NGC 1448 


31 


14(40) 


31 


15(54) 


31 


11(53) 


31 


13(27) 


2003ho 


ESQ 235-G58 


34 


10(41) 


34 


18(54) 


34 


13(53) 


34 


13(28) 


2003ip 


UGC 327 


33 


81(41) 


33 


75(54) 


33 


67(54) 


33 


76(28) 


2003iq 


NGC 772 


32 


50(40) 


32 


40(54) 


32 


34(53) 


32 


43(28) 


2004dj 


NGC 2403 


28 


06(43) 


28 


22(55) 


28 


19(54) 


28 


14(29) 


2004et 


NGC 6946 


28 


64(50) 


28 


21(58) 


28 


17(55) 


28 


37(31) 


2005CS 


NGC 5194 



















Note. — These values were calculated according to the approach exposed in 14.51 The mean value is an average of the values 
for each filter. 
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Table 4.5 EPM distances (|Jones et AL.IboOSi ) and SCM distances 



SN name 


dE96^ 


rl b 

"D05 


dscM 




(Mpc) 


(Mpc) 


(Mpc) 


1992ba 


16.4 (2.5) 


27.2 (6.5) 


18.5 (2.4) 


1999br 




39.5 (13.5) 


23.6 (3.4) 


1999eni 


9.3 (0.5) 


13.9 (1.4) 


10.0 (1.3) 


1999gi 


11.7 (0.8) 


17.4 (2.3) 


12.6 (0.6) 


2002gw 


37.4 (4.9) 


63.9 (17.0) 


48.1 (6.2) 


2003T 


87.8 (13.5) 


147.3 (35.7) 


108.1 (14.0) 


2003bl 




92.4 (14.2) 


65.2 (9.0) 


2003bn 


50.2 (7.0) 


87.2 (28.0) 


52.5 (6.8) 


2003hl 


17.7 (2.1) 


30.3 (6.3) 


25.6 (3.3) 


2003hn 


16.9 (2.2) 


26.3 (7.1) 


16.8 (2.1) 


2003iq 


36.0 (5.6) 


53.3 (17.1) 


30.6 (4.0) 



Eastman et al 



1 1991 



EPM distances from atmosphere models by 
^ EPM distances from atmosphere models by iDessart &: Hillierl ||2005| ) 



51 



2 - 



s 

O 1 .5 



T3 

I 



0.5 - 



-0.5 



"T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



d 



(D05) 



h-M- 



d 



(E96) 



H \ \ \ h 



H \ \ \ \ h 



H \ \ \ \ h 



-J-H 



H h 



J I I I I I I I I I I I I I I I I I I I I I I L 



1.2 



1.4 1.6 

log(dscM) 



1.8 



Figure 4.12 Direct comparison between SCM and EPM distances calculated by iJones et al.l (|2008f ). In 
magenta are shown the EPM distances computed from the atmosphere models of E96, while the green 
triangles are EPM distances obtained from the D05 models. The bottom panel shows the fractional differences 
between both techniques against \og{dscM)- In both panels the dotted lines trace the systematic shifts of 
the EPM distances with respect to the SCM distances, ~ 40% and ~ 12% using DOS and E96 atmosphere 
models, respectively. 
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Figure 4.13 Similar to Figure 14.121 after bringing the EPM distances to the SCM scale. The removal of the 
systematic differences leads to random differences between EPM and SCM distances. For the E96 case the 
discrepancies have a spread of 16%, while for the D05 case the scatter is 13%. 
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Chapter 5 



Discussion 



5.1 Variations of the extinction law 



The reddening law that minimizes the dispersion in the BVI Hubble di agrams (Ry = 1.4ifc0 .1) turns 



Cardelli at al 



19891 ). To our 



out to be very different than the standard Galactic extinction law {Ry = 3.1 
knowledge this is the first study of the interstellar extinction in external galaxies based on S Ne II-P. On the 
other hand, there have been several studies on this subject using Type la SNe. Most recently 



Folatelli et al 



(|2008l ) solved for Ry in a similar manner to our approach, i.e., by minimizing the di spersion in the Hubble 



diag ram and obt a ined a value of Ry 
and 



Reindl et al 



1.5, remarkably consistent with our value. 



Altavilla et al 



( 2004 ) 



( 20051 ) applied a similar procedure by minimizing the dispersion in the Mb — Ami^^B) 



relation for SNe la and found Rb ranging between 3.5 and 3 .7, significantly s maller than the standard 



value of 4.3. An even smaller value of Rb = 1-7 was found by ICapaccioli et al 



(Elias-Rosa et al 



( 1990 ) based on the same 



kind of ana. l ysis fo r SNe la. Further evidence for low Rx values were rep orted by 
Knot) et al. ( 2003)). Stu dies of indivi dual Type la SNe, such as SN 1999cl (iKrisciunas et al 



Phillips et all (|l999f ) and 



20071) . SN 2003cg 



20061 ) ■ SN 2002CV (|Ehas-Rosa et al.l 120081) . and SN 2006X (|Wang et al.l 120081 ) also yielded 
low values around Ry ~ 1.5. Other studies of dust reddening based o n diverse methods were performed to 



nearby galaxies obtaining val ues for Rb ranging between 2.4 and 4.3 (|Rifatta 



1992 



Brosch fc Loingei 



1990; 



Delia Valle fc Fanagia 



19911 ). The ra tio of total to selective ab sorption varies significantly between Ry 



2.6 — 5.5 even within our own Galaxy (jClavton fc Cardellilll988l ). A lower value of Ry could be due to dust 
grains smaller than those in our Galaxy, since for a given va lue of Ay the E{B — V) reddening decreases if the 
size of the grains grows. On the other hand. 



Wang 



(|2005f ) suggests the idea that scattering by dust clouds 
located in the circumstellar medium of the SN tends to reduce the effective R\ in the optical. This effect 
should be opposite in the ultraviolet, hence it could be further tested with photometry at these wavelenghts. 



5.2 Hq comparison with other methods 

The HDs constructed from the V~I color at day -30 as the extinction estimator (§ I4.3.ip give 
values of Hq between 70-73 km s^^ Mpc^^, which turn out to be very similar to those derived from the HDs 
where we leave Ry as a free parameter (70 ± 8 km s^^ Mpc^^). This shows that if we use V -I color based 
extinctions, the value of the Hubble constant is not be too sensitive to the adopted Ry. 
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Table 5.1 Hq values from the literature 



Method 


Ho 


Reference 




(km Mpc^i) 




SNe la 


71 ± 2 


Freedman et al. (20011 


SNe la 


71 ± 7 


Altavilla et al. f2004) 


SNe la 


73 ± 4 


Riess et al. f2005) 


SNe la 


62 ± 1 


Sandaee et al. (2006) 


SNe la 


72 ± 4 


Wang et al. (2006) 


SNe II-P (EPM) 


73 ± 6 


Schmidt et al. ( 1994) 


SNe II-P (EPM) 


52 ± 4 


Jones et al. (20081 


SNe II-P (EPMf 


92 ± 7 


Jones et al. (20081 


SNe II-P (SCM) 


55 ± 12 


Hamuv & Pinto (20021 


SNe II-P (SCM) 


72 ± 6 


Hamuv (20031 


SNe II-P (EPM) 


72 ± 9 


Freedman et al. (20011 


SNe II-P (SCM) 


70 ± 8 


this study 


TuUy-Fishcr 


71 ± 3 


Freedman et al. (20011 


TuUy-Fishcr 


59 ± 2 


Sandage et al. (2006) 


TRGB^ 


62 ± 2 


Sandage et al. (20061 


SBF'^ 


70 ± 5 


Freedman et al. (20011 


FP= 


82 ± 6 


Freedman et al. (20011 



TRGB = tip of the red giant branch l lSakai et al.ll2004h 

SBF = surface brightness fluctuations i Tonrv fc Schneider 
FP = "fundamental plane" 



1988 



* Using dilution factors of 
** Using dilution factors of 


Dessart & Hillier 


2005r 


Eastman et al. 


1 


1996) 



Table 15.11 summarizes several modern measurements of the expansion rate of the Universe, Hq, 
derived from different methods: SNe la, EPM and SCM for SNe II-P, the TuUy-Fisher relation, tip of the 
red giant branch (TRGB) distances, the surface brightness fluctuations (SBF) method, and the "fundamental 
plane" method for early type galaxies. According to Table 15.11 the value of the Hubble constant ranges 
between 52-82 km s^^ Mpc^^, and the most accepted value today is 70 km s^^ Mpc^^, which is in very 
good agreement with our SCM value of 70 ± 8 km s^^ Mpc^^. 

This study shows that the SCM can deliver accurate distances. The BVI Hubble diagrams yield 
scatter of 0.3 mag which implies a precision in individual distances ^ 15%. Part of this scatter could be 
due to the peculiar velocities of the SN host galaxies and the intrinsic precision of SCM could be even lower 
than this. In fact, when we perform a comparison for SNe in common between SCM and EPM, the distance 
differences range between 13-16% (after removing the systematic differencce among the SCM and EPM). This 
comparison is independent of the SN host-galaxy redshifts and implies that the internal precision in any of 
these two techniques must be lower than 13-16%, since these differences comprise the combined uncertainties 
of both methods. This is an encouraging result, since it implies that both EPM and SCM can produce high 
precision relative distances, thus offering a new route to cosmological parameters. 
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5.3 Final remarks 



We note that the dispersion in the HDs increases with wavelength. This seems contrary to expec- 
tations given that 1) the extinction effects decrease with wavelength, and 2) metallicity affects the B-band. 
Perhaps the luminosity is not only a function of velocity but also of metallicity, i.e., My — f{v, [Fe/H]). If so, 
the velocity term in eqs. l4.5H4.6l might be removing metallicity effects, more efficientlyat shorter wavelengths. 
This could be tested with [Fe/H] measurements of the Type II-P SNe. We plan to address this issue in the 
near future. Given the evidence we have, we can only claim a dependence of L with v. The obvious physical 
explanation is that the internal energy is correlated with the kinetic energy. The implication is that the ratio 
between the internal energy and kinetic energy is approximately independent of the explosion energy. 
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Chapter 6 



Conclusions 



We established a library of 196 SN II optical spectra and developed a code which allows us to correct 
the observed photometry for Galactic extinction, ii'-ternis and host-galaxy extinction. We applied our code 
to the BVRI photometry of 37 SNe II-P (§ 13. ip . We developed fitting procedures to the light curves, color 
curves and velocity curves which allow us to precisely determine the transition time between the plateau and 
the tail phases. The use of this parameter as the time origin permited us to line up the SNe to a common 
phase. The additional benefit of these fits is the interpolation of magnitudes, colors and velocities over a 
wide range of epochs. The methodology explained above yields the following conclusions: 

1. The compariso n between our color-based dereddening technique and the spectroscopic reddenings of 



Dessarti (|2008f) is satisfactory within the errors of both techniques. This is particularly encouraging 



since our method uses late-time photometric information, while the other method uses early-time spec- 
troscopic data, i.e. completely independent information. 

2. Using our new sample of SNe we recover the luminosity-velocity trend {LEV relation) previously 
reported bv lHamuv &: Pintd (j2002l ). 



3. We construct BVI HDs using two sets of host-galaxy reddenings. We demonstrate d that the (V 



Dessart 



/)-based extinctions do a much better job than the spectroscopic determinations by 
reaching dispersions of ~ 0.4 mag in the Hubble diagrams. This scatter is somewhat higher than that 
found previously bv iHamuvl (j2003l ) of 0.35 mag. A much smaller dispersion of 0.3 mag was achieved 
when we used V -I colors to estimate reddening and allowed Rv to vary. We obtain Ry — 1.4 ± 0.1, 
much smaller than the Galactic value of 3.1. The low value of Ry can be explained by a different nature 
of the dust grains in host-galaxies along the line of sight to Type II-P SNe. 

4. We derive a Hubble constant of 70 ±8 km s~^ Mpc^^ from BVI photometry, calibrating our HDs with 
Cepheid distances to SN 1999em and SN 2004dj, which agrees very well with the value obtained by the 
HST Key Project (jpreedman et allboOll ). 



Finally, we cal culate the d i stance moduli to our SN sample, and make a comparison against EPM 
distances from 



Jones et al 



()2008l) . The 11 SNe in common show a systematic difference in distance 
between EPM and SCM, depending on the atmosphere models employed by EPM. Correcting for these 
shifts we bring the EPM distances to the SCM distance scale, from which we measure a dispersion of 
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13-16%. This spread reflects the combined internal precision of EPM and SCM. Therefore the internal 
precision in any of these two techniques must be < 13-16%. 

This analysis reconfirms the usefulness of SNe II-P as cosmological probes, providing strong encour- 
agement to future high- 2; studies. We found that one can determine relative distances from SNe II-P with a 
precision of 15% or better. This uncertainty could be further reduced by including more SNe in the Hubble 
flow. In its current form the SCM requires both photometric and spectroscopic data. Since the latter are 
expensive to obtain (especially at high-^) it would be desirable to look for a photometric observable as a 
luminosity indicator instead of the expansion velocities. 
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Appendix A 



The Computation of Synthetic 
Magnitudes 



The implementation of SCM requires the implementation of AKA corrections to the observed SN 
magnitudes (§ I3.ip . This process involves synthesizing broadband magnitudes from the library of SNe ITP 
spectra. It is crucial, therefore, to place the synthetic magnitudes on the same photometric system employed 
in the observations of the SN. 

Since the SN magnitudes are measured with photon detectors, a synthetic magnitude is the convo- 
lution of the observed photon number distribution, Nx , with the filter transmission function 5'(A), i.e.. 



mag 



2.5 logiQ / Nx S{X) dX + ZP, 



(A.l) 



where ZP is the zero point for the magnitude scale and A is the wavelength in the observer's frame. 

For an adequate use of eauation lA.il S(A) must include the transparency of the Earth's atmosphere, 
the filter transmission, and the detector q uantum efficiency (QE). For BVRI I adopt the filter functions -Bgo, 
Vgo, -R90, -^90 pubhshed bv iBesselll (|l990l ). However, sin ce these curves are meant for use with energy and 



not photon distributions (see Appendix in 



Bessell 



19831 ) , I must divide them by A before employing them in 



equation lA.il Also, since these filters do not include the atmospheric telluric lines, I add these features to the 
R and / filters (in B and V there are no prominent telluric features) using my own atmospheric transmission 
spectrum. Figure lA. II shows the resulting curves. 

The ZP in cq. lA.ll must be determined by forcing the synthetic magnitude of a sta r to m atch its 
observed magnitude. I use the spectrophotomctric calibration of Vega published by HavesI (|l985r ) in the 



Johnson et al 



(jl966l ) , from which I solve 



range 3300-10405 A and the V magnitude of 0.03 mag measured by 
for the ZP in the t^-band. In principle, I can use the same procedure for BRI, but Vega's photometry 
in these bands is not very reliable as it was obtained in the old Johnson standard system . To avoid these 



problems I employ ten stars with ex cellent spectrophotometry (jHamuv et al.lll992. 19941 ) and photometry 



in the modern Kron-Cousins system (jCousins 



1971 



198C 



19841) . Before using these standards I remove the 



telluric lines from the spectra since the filter functions already include these features. With this approach 
I obtain an average and more reliable zero point for the synthetic magnitude scale with rms uncertainty of 
~0.01 mag. With these ZPs I find that the synthetic magnitudes of Vega are brighter than the observed 
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Figure A.l BVRI filters functions of lBesselll ([19901) meant for use with energy distributions (dotted curves). 
With solid lines are shown the curves modified for use with photon distributions, to which I added the telluric 
lines. 
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Figure A. 2 Adopte d spectrophotometric calibration for Vega. In the optical (A < 10,500 A) the calibration 



is from lHavesI (|1985I ). and at longer wavelengths I adopted the Kurucz spectrum with parameters Te/y=9,400 
K, log .g=3.9, [Fe/H]=-0.5, V„,^croturb=0. 



61 



magnitudes (jjohnson et al.lll966r ) by 0.016 mag in B, 0.025 in R, and 0.023 in / (see Table \K7l^ . which is 
not so surprising consi dering that t his comparison requires transforming the Johnson RI magnitudes to the 
Kron-Cousins system (jTavloiill986l ). Figure E2] shows the adopted spectrophotometric cahbration for Vega. 



Table A.l Photometric Zero points and Synthetic Magnitudes for Vega 





B 


V 


R 


/ 


Zero point 


35.287 


34.855 


35.060 


34.563 


Vega 


0.014 


0.030 


0.042 


0.052 



Table ETT] summarizes the zero points computed with eq. lA.ll and the corresponding magnitudes for 
Vega in such system. For the proper use of these ZPs it is necessary to express N\ in sec~^ cm~^ cm^^ and 
A in A. From the ten secondary standards 1 estimate that the uncertainty in the zero points is ^ 0.01 mag 
in BVRI. 
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Appendix B 

Finding J'{v) using the DSM in 
Multidimensions 



Given the S-parameter function we have to deal with (refer to J- in eq. I3.2.1|) . we had to ex- 
plore different method s of multidimensional minimization. The Downhill Simplex Method (DSm]^, due to 



Nelder fc MeadI (|1965[ ). does not use one-dimensional minimization as a part of their computational strategy. 
Derivatives are not required, only function evaluations. Its funcionality is based on the simplex, a geometrical 
figure consisting in N dimensions, of iV+1 vertices, all their interconecting line segments, and polygonal faces. 
For example, in 2-dimensional space a simplex is a triangle. The starting guess should not be just a single 
point, but + 1, defining an initial simplex. The DSM now takes a series of steps (shown in Figure [B.l|) . 
most steps just moving the vertice of the simplex where the function is largest through the opposite face of 
the simplex to a vertice with a lower value of the function. The routine amoeba, called after the unicellular 
organism, tries to "swallow" the minimum being intended to be descriptive of this kind of behavior. The 
hungry amoeba is the simplex which is "fed" with minima. The purpose after a few step is the reduction 
of the TV-dimensional volume of the simplex as shown in Fig. IB. II for the 3D case. Hence the minimum gets 
enclosed and finally found with an accuracy given by the size of the simplex. 

The input parameters for the Fortran version of this routine are the following 

amoeba(isimplex[mp,np] ,y[mp] ,mp,np,ndim,ftol, funk, iter) 

where 

isimplex[nip,np] : the matrix containing the coordinates of the initial simplex, np+1 vertices defined by a 
np-dimensional vector 

y[mp] : vector composed by the value function evaluated in each of the np+1 initial vertices (rows) of 
isimplex 

mp = np + 1 

np = ndim 



^Consult for further details Numerical Recipes in Fortran (chapter 10.4: lPress et al1ll992t) as the main bibliographic reference 
of this appendix. 
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Figure B.l Steps of the amoeba algorithm trying to "swallow" the minimum using the DSM. This example 
draws the simplex in 3-D (tetrahedron). The possible outcomes, that move away from the highest point, are 
(a) the vertice is reflected through its opposite face, (b) a reflection plus a 1-D expansion, and (c) a 1-D 
contraction. A more unusual step performs (d) a contraction along all dimensions towards the lowest vertice. 
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ndim ; number of dimensions or variables of the function (parameters to be fitted) 
f tol : the fractional convergence tolerance to be achieved in the value of the function 
funk ; function of the form funk(x) where x [mp] is a np+l-dimensional vector 
iter : total number of function evaluations after the convergence 

The output consists of the isimplex matrix and the y vector overwritten by np+1 new points all 
within a tolerance f tol of the minimum function value, whereas the number of function evaluations is saved 
in iter. Thus, the result is a new and smaller simplex containing the minimum. This procedure guarantees 
a succesful search of at least a local minimum value of the function. 

In our specific case we want to minimize x^i"^) ^ function of the parameter vector v, which is 
written as 

xHv) - E ( "^°gH--^(f'^H) ') ' (B.i) 



where T is the function defined in eg. 13. 81 whose role is to model the SNe II-P hght curves, [magli] ± CTmagH) 
is the observed magnitude of the SN with its associated uncertainty for a given Julian date t[i], all together 
completing a total of nobs measurements that, basically, draw the observed SN light curve. 

Although simple and user-friendly, the DSM finds problems working on an 8-dimensional space but 
not as most algorithm would do. In most cases the initial simplex has to be accurately defined according 
to the desired results. In order to make a semi-automatic search for the initial simplex, we introduce new 
parameters in the minimization routine. From the photometric observations we extract the date of the first 
data point ti, and using a simple algorithm we calculate the approximate end of the plateau tend (refer 
to § 13.2. ip . By means of these calculations we constrain the J-'-parameters having temporal components. 
In order to contract the magnitude component of the J-"-parameters we extract from the photometry the 
brightest and the dimmest magnitude, m{,ri and rridim respectively. For example, below is shown the initial 
simplex we feed into amoeba to fit the V light curve of SN 1999em, 



isimplex = 



ao 


tpT 


Wo 


Po 


mo 


P 


Q 


R 


( ^ 


tend 5 


4.5 


0.08 


nibri 


-1.2 


60 


t^ + 9 \ 


1.6 


tend ~t" 25 


45 


0.08 




-1.2 


60 


tend — 10 


1.6 


tend 5 


8 


0.08 




-1.2 


60 


tend — 10 


1.6 


tend 5 


4.5 


0.012 


rribri 


-1.2 


60 


tend — 10 


1.6 


tend 5 


4.5 


0.08 


25 


-1.2 


60 


tend — 10 


1.6 


tend 5 


4.5 


0.08 




1.7 


60 


tend — 10 


1.6 


tend 5 


4.5 


0.08 


rribri 


-1.2 


21 


tend — 10 


1.6 


tend 5 


4.5 


0.08 


rribri 


-1.2 


60 


tend — 10 


I 4 


tend ~t- 25 


8 


0.012 


25 


1.7 


60 


t^ + 9 J 



(B.2) 



where the rows are the 9 initial vertices of the simplex in the 8-parameter coordinate system. This is one 
structure of the many one can choose to build the isimplex that gives us the best convergence. Each filter 
requires slightly different parameter ranges. By averaging the vertices of the final converging simplex, we 
obtain a unique parameter vector with its corresponding uncertainty vector calculated from the dispersions 
of the 9 vertices around the mean. The resulting fits for SN 1999em are shown in the upper-right panel of 
Fig. 13.41 and the corresponding parameters are summarized in Table IbTT] 
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Table B.l J^-parameters for the V light curve of SN 1999em 



ao 


1.978 ±0.015 


mag 


tpT 


JD 2451590.1 ±1.1 




Wo 


4.52 ±0.85 


day 


Po 


9.0 ±0.4 


10~^ mag/day 


mo 


16.297±0.014 


mag 


P 


0.403 ±0.019 


mag 


Q 


JD 2451475.3 ±3.7 




R 


64.6 ±5.7 


day 
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